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Abstract. We study an inverse problem for a non-compact Riemannian man- 
ifold whose ends have the following properties : On each end, the Riemannian 
metric is assumed to be a short-range perturbation of the metric of the form 
(dy)'^ + h{x, dx), h{x, dx) being the metric of some compact manifold of codi- 
mension 1. Moreover one end is exactly cylindrical, i.e. the metric is equal to 
(dy)-^ +h{x, dx). Given two such manifolds having the same scattering matrix 
on that exactly cylindrical end for all energy, we show that these two manifolds 
are isometric. 

1. Introduction 

The aim of this paper is to study spectral properties and related inverse problems 
for a connected, non-compact Riemannian manifold of dimension n > 2 with or 
without boundary. We assume that is split into -I- 1 parts 

(1.1) r2 = /cur2iU---ui7Ar, 

where /C is a bounded open set, and fi^, called an end of 11, is diffeomorphic to 
Mi X (0,oo), Mi being a compact manifold of dimension n — 1. (See the figure 1.) 
Denoting the local coordinates on Mi by x, we assume that Mi is equipped with a 
Riemannian metric hi{x, dx) = Y^p qLi hi,pq{x)dxPdx'^ . Letting y be the coordinate 
on (0, oo), we denote the local coordinates on Hi hy X — {x,y). We assume that 
the Riemannian metric G on fi, which is denoted by Gi — ^^j^ gi,pq{X)dXP dX'' 
on fli, has the following property 

(1.2) |9^(g,,p,(X)-/i,,p,(x))| <C„(l + y)-i-^«, Va, 

where hi^pn{x) = hi^np{x) = Oifl < p < n — 1 and hi_nn{x) = 1, and Cq is a 
constant. The metric hi{x,dx) on Mi is allowed to be different for different ends. 
We shall assume either il has no boundary or each Mi, consequently itself, has a 
boundary. In the latter case, we impose Dirichlet or Neumann boundary condition 
on dft. Let H = — Ac where Aq is the Laplace-Beltrami operator associated with 
the metric G. One can then define a scattering operator S{X) = (S'ij(A)), which is 
a bounded operator on L^(Mi) © • • • © LF'{Mn), where A g {Eq^oo) \ £{H) is the 
energy parameter, Eq ~ mi (Tess{H), and £{H) is the set of exceptional points to 
be defined in (|3.34p . Our goal is the following. 



Date: May 11, 2009. 



2 



HIROSHI ISOZAKI, YAROSLAV KURYLEV, AND MATTI LASSAS 



^2 




Figure 1. Manifold has ends fi^, j = 1, 2, . . . , N. 

Theorem 1.1. Suppose we are given two manifolds r = 1,2, of the form 

il.l]) having Nr ends, i — 1, - ■ ■ , iV^, equipped with the metric G^^^ satisfying 

the assumption U.^^) . Assume that fi^^^ = fl^^^ and 

ri-l 

(1.3) G^^'^ =G^^'' = {dyf + hi{x,dx), hi{x,dx) ^ ^ hi,jk{x)dx^ dx^ 

onn[^'' =n^^\ moreover S[\\X) = S[f {X) for all X G {E' ,oo)\{£^^^ US'^^^), where 
f (""^ is the set of exceptional points for H^'^\ and E' = max (£;^^\£;^^^). ThenQ^^^ 
and are isometric as Riemannian manifolds with metrics G'^^\ G^^^. 

This means that if we observe waves coming in and going out of one end fii, 
which is assumed to be non-perturbed, we can identify the whole manifold Q. Note 
that in Theorem 11.11 neither the number of ends of each fJ^'') nor the metric on 

(r) 

the manifold are assumed to be known a-priori. The key idea of the proof is 

to introduce generalized eigenfunctions of the Laplace-Beltrami operator which are 
exponentially growing at infinity, and define the associated non-physical scattering 
amplitude. The crucial fact is that this non-physical scattering amplitude is the 
analytic continuation of the physical scattering amplitude. Then the physical scat- 
tering amplitude determines the non-physical scattering amplitude, which further 
determines the Neumann-Dirichlet map of the interior domain. By the boundary 
control method (see [31 [SI [71 [311 ISZ] ) , one can determine the metric inside. 

In this paper, we exclusively deal with the Neumann boundary condition. The 
other cases are treated similarly and in fact more easily. The forward problem, 
i.e., the existence and completeness of wave operators, the eigenfunction expansion 
and the representation of S-matrices are well-known for short-range perturbations 
(see e.g. [l^, [23], [H], [12], [SU], [M], [H], see also [13]). The new issue we have 
to discuss in this paper is the difference of conormal derivatives on the boundary 
associated with unperturbed and perturbed metrics. Therefore, focusing on this 
point, we only explain the outline of the proof of the forward problem under the 
assumption (|1.2p following the approach in 30J, where spectral theory and inverse 
problems on hyperbolic spaces are developed in an elementary way. 
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In the cylindrical ends, the physical generalized eigenfunction of the Laplace- 
Beltrami operator admits the analytic continuation with respect to the energy pa- 
rameter, and this analytically continued eigenfunction is exponentially growing as 
y — » oo. This sort of non-physical exponentially growing generalized eigenfunction 
was first introduced by Faddeev to develop the multi-dimensional Gel'fand-Levitan 
theory (|19p. The exponentially growing solutions of Schrodinger equation was 
rediscovered in 1980's, where the new feature is its utility in the inverse scattering 
from a fixed energy as well as in the inverse boundary value problem (see [IS] , [17] , 
[44j . [55] and the review |t20j). The interesting fact is that this apparently mys- 
terious exponentially growing generalized eigenfunctions appear naturally in the 
cylindrical domain. Using these exponentially growing eigenfunctions, it is possible 
to obtain, from the scattering operator S'ii(A), the Gel'fand spectral data on a part 
of the boundary T = Mi x {1} of the non-compact manifold fl^ = ^l\ (Mi x (1, oo)). 
The Gel'fand boundary data for this case is the family of the Neumann-Dirichlet 
map, A.{z), A.{z)f = u\ , where u is the solution to the boundary value problem 



To solve this problem, we use the boundary control (BC) method (see [^ for the 
pioneering work and [7], [32] for the detailed exposition). We note that typically 
the BC method deals with inverse problem on compact manifolds. The case of 
non-compact manifold considered here requires substantial modifications into the 
method, since the spectrum is no more discrete and it is also impossible to use 
eigenfunctions as coordinate functions. A short description of the BC-method for 
non-compact manifolds was given in |34j. Here we provide detailed constructions 
for the considered case of a manifold with asymptotically cylindrical ends. 

The structure of this paper is as follows. Sections 2, 3, 4 are devoted to a detailed 
analysis of scattering on manifolds with asymptotically cylindrical ends. After some 
preliminary estimates for the case of a half-cylinder with a product metric in §2, we 
discuss the spectral properties of the Laplacin in in §3. Using these properties, 
we develop the scattering theory for such manifolds in §4. The remaining part of 
the paper. Sections 5, 6 are devoted to the inverse scattering. In §4, we show that 
S'ii(A) determines the Neumann-Dirichelt map A(z). An important step, which 
at the moment requires the product structure of the metric on Mi x (0, oo), is 
the recovery, from physical scattering matrix S'ii(A), the non-physical scattering 
amplitude. At last, §6 is devoted to the development of the B-method for non- 
compact manifolds. For the convenience of the reader, interested predominantly in 
the inversion methods, we make this section independent of the previous ones. 

Our manifold O is a mathematical model of compound waveguides, e.g. settings 
of optical and electric cables, oil, gas and water pipelines, etc, which are the most 
typical geometric constructions encountered in the every-day life. As for the inverse 



(-Agi - z)u = in n^, 
d^u = on dnndn\ 
dyU = dyU / on r. 
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problem, many works have been devoted so far to the distribution of resonances for 
the waveguides ([5], [H], [5], [Tl], [TS], [1], [IHI)- Identification or reconstruction of 
the domain or the medium for grating, layers or waveguides are studied by [12], [29] . 
[46j . [18] . In particular, a similar inverse problem for waveguides was considered by 
Eskin-Ralston-Yamamoto il8j when is a slab, (0, B) x K, with the variable sound 
speed c{x,y), where c{x,y) = c{x) for large \y\. The present paper deals with the 
forward and inverse scattering problems for waveguide in a full generality. 

The notation in this paper is standard. For a self-adjoint operator A, u{A), 
crp{A) and a^ssiA) mean its spectrum, point spectrum and essential spectrum, 
respectively. For two Banach spaces 7^1,7^2, B(7ii;7i2) means the space of all 
bounded operators from Hi to Ti.2- For an operator A on a Hilbert space H, D{A) 
denotes its domain of definition. For a Riemannian manifold A4, H"^{A4) denotes 
the usual Sobolev space of order to on A^. For a domain D and a Hilbert pace TC, 
L^{D;TC;dfi) means the space of 7i-valued L'^-functions on D with respect to the 
measure dfi. If 7i = C, we omit it. For a differentiable manifold M and p e M, 
Tp{M) denotes the tangent space of M at p. A simplified version of our results is 
given in [?T| . 



As a preliminary, let us begin with proving some a-priori estimates for the op- 
erator —dy — Afi on r^o = M X M_|_ with Neumann boundary condition, where 
y € R+ = (0,00), M is a compact Riemannian manifold, and A/j is the Laplace- 
Beltrami operator associated with metric h{x, dx) equipped on M. 

2.1. Besov type spaces on cylinder. We define an abstract Besov type space, 
which was introduced by Hormander 1 in the case of M". Let M be the above 
mentioned compact manifold, and (, )m, \\ ■ \\l^{m) be inner product and norm of 
L'^{M), respectively. We define intervals /„ by 



Let B be the Banach space of L^(Af)-valued functions on (0,00) equipped with 
norm 



2. A-PRIORI ESTIMATES IN HALF-CYLINDERS 





Its dual space is the set of 



(Af)-valued functions u{y) satisfying 
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It is easy to see that there exists a constant C > such that 



sup 2-/2 (^^ ll^(y)lli.(M)^2/)'^'< hup;^^''ll«(2/)lli^(M)^ 



1/2 



<C sup 2-/2 ^ \\v{y)\\h^M)dy 



1/2 



Therefore, we identify B* with the space equipped with norm 

1/2 

/ 1 1'- o , \ 

I u Lb* 



We also use the following weighted space and weighted Sobolev space: For s e R, 

/•CX) 

L''' ^ f ^ ml = / (1 + yfimWhmdy < oo, 
Jo 

iJ™'* 9 U <^ II'uIIh'".- = 11(1 + yyu\\H'n{Mx{0,oo)) < oo. 

In the following, || • || means || • ||o and (•, •) denotes the inner product of L'^(M x 
It often denotes the coupling of two functions / G L^'" and g G L^'^" or f E B 
and g € B* . The following inclusion relations can be shown easily, and the proof is 
omitted. 

Lemma 2.1. For s > 1/2, we have 

^2„, ^ g ^ ^2,1/2 ^ ^2 ^ ^2,-1/2 ^ ^* ^ 

We often make use of the following lemma, whose proof is also elementary and 
omitted. 

Lemma 2.2. Suppose u E B* . Then 

(2.1) }^^^j^''\Hy)\\h^M)dy = o, 

if and only if 

(2.2) p{^)\\<y^\\UM)dy = o, '^p^^o°m'^))- 

2.2. A- priori estimates. Let us consider the following equation in i^o = M x ]R_|_: 

r (-92 -Ah-z)u = f in Oo 
1 d^u = on dQo, 



(2.3) 



z being a complex parameter, and di, conormal differentiation on the boundary. In 
the following, we often denote by the norm of derivatives of |a|-th order of 

u without mentioning local coordinates. 

Lemma 2.3. Let z G C fee given. Then : 
(1) If u, f € 1/2. s j^Qj. gQfi2(, s (eR, we have 

^ R"a>|u<c(|HU + 11/11,). 

|a|+i<2 
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(2) If u, f B* , then we have 

\\dM\B' + \\dyU\\B' <C{\\u\\B^ + \\f\\B'). 

Proof. We shall prove (2). Pick xiv) & C^(M) such that xiv) = 1 ilvl < 1), 
x{y) — {\y\ > 2) and put Xniv) ~ xiv/R)- We take the inner product in L^(r2o) 
of (|2.3p and x\{y)'U- We then have 

2 y 

WXRdyuW^ + {xRdyU, -^X'{^)u) + WxndxuW'^ ~ zWxruW^ = {f,XRu), 
which implies 

WXRdyuf + WxRAuf < C + WxRuf + |lxfl/f ) • 

Then we have for i? > 1 

rR nR 

\\dyu\\l2andy + / \\dxu\\l2ai)dy 
Jo 

< C (^J^'' \\u\\h^M)dy + f^"^ ll/lli.(M)rf2/^ ■ 

Dividing by R and taking the supremum with respect to R, we obtain (2). 

Let us prove (1). The 1st order derivatives are dealt with in the same way as 
above. We put w = (1 + yYu. Then v satisfies {—dy — A/j — z)v — g, where 
g e -L^(il). By the a-priori estimates for elliptic operators, we have v G H'^(fl), 
which proves (1). □ 

Let Ai < A2 < • • • 00 be the eigenvalues of —A^, and P„ the associated 
eigenprojection. Then 

00 

\/z + Aji ^ ^ \/z- Xn Pn, 
n=l 

where for ( = re*", (r > 0,0 < 6 < 27r), we define = V^e'^/^. 

Our next aim is to derive some a-priori estimates for solutions to the equation 
We use the method of integration by parts due to Eidus (LI6J). We put 



P{z) = ^/z + Ah, 

D±{z) = dyTiP{z). 
Then the equation (|2.3p is rewritten as 

(2.2) dyD±{z)u = TiP{z)D±{z)u~ f. 

Lemma 2.4. Let (p{y) G C°°(M) be such that (p{y) > 0. For a solution u of the 
equation h2.S\) . we put w — D-^-{z)u. Then iflmz > we have for any < a < b < 
00 

y=b 



f'{y)\\w{y)\\L2(M)dy <2 ^iy)\if,w)L2(M)\dy + '^||w||l2(m) 
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Proof. Since w satisfies dyW = —iP{z)w — /, we have 

Lp{y){dyW,w)L2(M)dy ^ -i j ^p{y){P{z)w,w)L'^(M)dy- j 'f{y){f,w)L^(M)dy. 

•J a J a 

Integrating by parts and taking the real part, we have 

1 b l-b 

v>\\w\\l2^M) - / ^'iy)\\wiy)\\h(M)dy 

ip{y){lmP{z)w,w)L2(^M)dy - 2Re / (p{y){f,w)L2i^M)dy- 

J a 

Taking notice of ImP(z) > for Imz > 0, we get the lemma. □ 

Let C+ = {z e C;Imz > 0}. 
Lemma 2.5. Let w be as in Lemma 2.4 and suppose that 

1 

(2.3) IT^rJ^ \Hy)\\hiM)dy = o. 

Then there exists a constant C > independent o/ z G C+ such that 

\Hy)\\l2^M)<C\\f\\B\\w\\B', VyeM. 
Proof. Taking (p{y) = 1 in Lemma 2.4, we have 

.b 

lk(a)lli2(M) < l|w^(&)lli2(M) + 2 / \if,w)L^M)\dy 

J a 

< \\wmhiM)+C\\f\\B\\w\\B'. 

The assumption of the lemma implies liminf ||i(;(6)||^2/jvf) — 0, which proves the 
lemma. □ 

Corollary 2.6. Under the assumption of Lemma 2.5, there exists a constant C > 
such that 

\\w\\B*<C\\f\\B, yzec+. 

Proof. Lemma 2.5 implies that 

1 /"^ 

Iklle- = sup - / \\w{y)\\l2andy < C||/||b||w||b., 

R>1 ^ Jo 

which proves this corollary. □ 

Theorem 2.7. For a small 6 > 0, let 

Js = {z e C+ ; dist {Rez,a{-Ah)) > 5}. 

Let u be a solution to i2.3\) such that w = D^{z)u satisfies i2.3\} . Then there exists 
a constant C > such that 

\\u\\B^<C\\f\\B 

holds uniformly for z G . 
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Proof. Let A{z) = ReF(z) (P(z) + P{z)*)/2. By the equation we have 

In view of the formula 

-i{P{z)w,u)L2(M) = -2i{A{z)w,u)L2(M) +i{P{zyw,u)L^M) 
= -2i{w,A{z)u)L2(M) + i{w,P{z)u)L^M), 

we then have 

dy{w,u)L2(M) = -2«(w, A(z)u)l2(m) - {f,u)L2(M) + \\M\h{M)- 

Using w — dyU — iP{z)u, we compute 

2i{w,A{z)u)L2(^M) = '^i{dyU,A{z)u) L2f^^i) + \\P{z)u\\\2(m) + {P(zfu,u)L2^M)- 
Summing up, we have arrived at 

dy{w,u)L2^M) = -'^i{dyU,A{z)u)L2(M) - \\P{z)u\W2^^m) 

- ((z + A,Jm,m)l2(m) - (/, m)l2(m) + ||w||^2(M)- 
Taking the imaginary part and integrating in y, we have 



Im 



(w,u)l2(m) 



y=b 



y=a 



2Re / {dyU, A{z)u)l2(^m) 



Imz / ||u||^2(M)Of2/ - Im / (/, w)L2(Af)dy. 
<y a <y a 



Since A{z) is self-adjoint, we have by integration by parts 



2Re 



{dyU,A{z)u)L2^M)dy = {A{z)u,u)l2(^m) 



y=b 



Using Im z > 0, we obtain 



(2.4) 



Im 



{w,u)l2(m) 



y=b 



y=a 



1 V=b 

(A(z)u, M)i2(M) < C||/||b||u||b-, 

J y=a 



where C is independent of z G <C+. We renumber the eigenvalues of — A/i in the 
increasing order /xi < /i2 < • • ■ without counting multiplicities and put //q = — oo, 
i.e. {A„ ; n = 1, 2, • • • } and {/i„ \ n — 1, 2 • • • } are the same as subsets of R. For a 
sufficiently small (5 > 0, we put 

Jn,s = {z e C+ ; + ^ < Re z < ^„ - (5}. 

Assume z S Jn,5 and split uasu = M<+u>, where 

U< = ^ PjU, M> = ^ PjU, 

Recall that Pj is the eigenprojection associated with Xj . We also define w< , w> , /< , 
/> similarly. Note that w;< = D+(2;)u<. Let us remark that (|2.3p and therefore 
(|2.4p hold with w,u,f replaced by w^,u^, and Wy,u^,f^, respectively. For 
eigenvalues \j < /in-i, we have Re z ~ Xj > \/6. Therefore 

(2.5) (A(z)M<,u<)i2(M) > V^|lu<|li2(M)- 
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Since dyu{0) = 0, we have w<(0) = — «P(z)u<(0). Therefore 

-Im (u;< (0), M<(0))i2(M) = Re (P(z)m< (0), u<(0))l2(m) 

= (A(2;)u<(0),m<(0))l2(m). 
Letting a = 0, b = t in (12. 4p . we then have 

Im(u;<(i),u<(t))i2(M) + (^(z)w<(i),u<W)L2(M) < . 
Using (|2.5p . we have 

b\\u<\\b')- 

Using Corollary 2.6, we then have 

^ £ \\u<iy)\\h^M)dy < C (|1/<||2 + WUWbMb') , 

which implies 

(2.6) \\u<\\b' <C\\f<\\B. 

On the other hand, if Xj > /x„, we have Re (Xj — z) > S. Therefore 

(2.7) {~dl - A;,. - z)u> = i-dl + ilmz)uy = />, 

where is a uniformly, with respect to z, strictly positive operator on L^(AI). 
Hence, we have 

(2.8) II«>IIl^ <qi/>l|L^, 

which by Lemma 2.1 implies 

(2.9) \\u>\\b' <C\\f>\\B. 

The above two inequalities (|2.6p and (|2.9p prove the theorem. □ 

3. Manifolds with cylindrical ends 

3.1. Resolvent equation. We return to the manifold fl = /C U fii U • . . U flN 
introduced in §1. Fix a point Pq £ IC arbitrarily, and let dist(P, Pq) be the geodesic 
distance with respect to the metric G from Pq to P. We put 

noiR) = {Pen; dist (P, Po) < R}, n^iR) = {P e O ; dist (P, Po) > R}. 

For P > large enough, take xo G C'o°(^) ^'^'^h that xo = 1 on fio(^)7 Xo = 
on Qoo{R +1). Define Xj ~ ^ ~ Xo on fij, Xj = on \ fij. Then {xjljlo ^ 
partition of unity on fi. 

Let Ag be the Laplacc-Beltrami operator for the metric G on fl endowed with 
Neumann boundary condition on dfl. The conormal differentiation with respect to 
G is denoted by d^. We put 

H = -Ag, R{z)^{H-z)-\ 

As in §1, we identify fij with Mj x (0,oo), and let hj{x,dx) be the metric on Mj. 
We compare G with the unperturbed metric G^^^ = (dy)^ + hj{x,dx) on Qj. Let 



10 



HIROSHI ISOZAKL YAROSLAV KURYLEV, AND MATTI LASSAS 



(3.1) 



A„(o) be the Laplace-Beltrami operator for with Neumann boundary condition 
on 9r2j . The associated conormal differentiation is denoted by d m . We put 

Our next concern is the difference between the boundary conditions for H and 
H^°\ We put for large i? > 

Lemma 3.1. There exists a real function w{x^y) G C°°{U,j) such that 
d„w{x,y) = on dQj{R) 
w{x,y) — y + 0{y^^^'^°) as y ^ oo. 

Proof. By the decay assumption (|1.2p . letting w{x,y) — y + w{x,y), we should 
have d,yW = —d,yy = 0{y~^~'^°) on dil,j{R). Extending the vector field v near the 
boundary and integrating along it, we get w — 0{y^^^'^"). □ 

For m > and s e M, we define the weighted Sobolev space on the boundary by 

^ e i7™'^(9%(i?)) ^ (1 + yY^j e H"\dQ,iR)). 

Lemma 3.2. There exists an operator of extension £j such that for m > 1/2 and 

ijj g H™(a%(i?)) 

> on 9%(i?), 
on f7\(%nf7oo(i?-l/2)), 



(3.2) d,£jij 



(3.5) 



(3.3) supp (Sjip) c ilj n Ooo(i? - 1). 
For m > 1/2 and s > 0, it satisfies 

(3.4) £j e B(i/"'''(9%(i?));iJ"+3/2-^(17j-)). 

Proof. Let M' = flj n flQo{R — 2). We smoothly modify the corner of M' , i.e. 
{P e n dQ ; dist(P, Pq) = ^ - 2}, and let M be the resulting manifold. Let vm 
be the unit outer normal to A4 . By solving the elliptic boundary value problem 

f {-AG + i)u = in M, 

\duMU = i^ on dM, 

we define Ejijj = x^m, where Xj G C'°°(rij) is such that Xj = 1 on ^jOrtooiR— 1/4), 
Xj = on f]\ (rjj nrJoo(i?- 1/2)). it then satisfies dSm, ([SSI)- The property ((O)) 
for s = follows from the standard estimate for the elliptic boundary value problem. 
Let < s < 1 + eo and take ip G H™'^{dM). For the solution u to the boundary 
value problem p.5p . we put m = {l + w{x,y))~'^u' and "0 = (1 + w;(x, y))^^?/'', where 
w{x,y) is constructed in Lemma 3.1. Then u' is a solution to the boundary value 
problem 

r(-AG + i' + K)?i' = in M, 
\di,j^u' = i}" on 
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where k > is sufficiently large, and L' is a 1st order differential operator with 
bounded coefficients, and ip" = tp' on dQj{R). Since the mapping ■(/;" — s- u' is 
bounded from H'^idM) to H''^+^/'^{M), we get dM]) with < s < 1 + eq- Repeat- 
ing this procedure, we can prove (|3.4p for all s > 0. □ 



For u e H^{Q.j) satisfying d (o)M = on dilj{R), we have 

(3.6) d^{xju)^w{x,yy^-^"BjU on dQjiR), 
where 

(3.7) = w{x,yy+'" {xjid. - 5^(o,) + (d^Xj)) 

is a 1st order differential operator on dilj{R) with bounded coefficients. We put 

(3.8) £j ^w{x,yr^-'">£j. 

Then by jS^]) and dSH), for u e i72(f]) satisfying d^mu = on dn^iR) the 

following formula holds 

(3.9) d^SjBjU^ d^ixju) on af7j(i?). 
Moreover 

(3.10) y^+'°£jBj e B(ij2(17);i/2(f]))nB(i7^/2(^).^3/2(^))_ 
Suppose u satisfies 

r (-Ag(o) ~ z)u^ f in rij, 
|(9^(o)U = on fijndn. 

Then by (|3.9[) . tij = XjU ^ ^jBjU satisfies 

f (-Ag-z)wj =Xj7 + V, (2)u in 



(3.11) 



=0 on 5r2, 



where 

(3.12) Fj(z) = [-Ag,Xj] + Xj(Ag(o) - Ag) + (Ag + z)£:,Bj. 

Lemma 3.3. Lei G C°°(il) &e smc/i i/iat Xj — \ on Vlj n iloo{R — 1) a?Jc? Xj 
outside flj n r2oo(-R — 2). T/ien /or z ^ R, i/ie following resolvent equations hold : 

(3.13) i?(z)x, = (x, - S,B, ~ R{z)V,{z))Rf\z)xj. 

(3.14) XjR{z) = XjJ-'Rf{z)Jj {xj - {SjBjY - V,(z)*i?(z)) , 

where Jj ~ (dot G/det G^*''')"'^''^, and the adjoint * is taken with respect to the inner 
product of L'^iyi) with volume element from the metric G. Moreover R^^ {z)Jj{£jBj)* 
ana Rf\z)JjVj{z)*R{z) are compact on L^(f2). 
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Proof. Let u = Rj {z)xjf for z ^ R. Then checking the boundary condition 
by (EH), we have Vj = X]Rf\z)X]f - £-jB,Rf\z)xjf e D{H), and by (EU) 
{H — z)vj = XjXjf + Vj{z)u = Xjf + Vj{z)u, which imphes (|3.13p . 

By extending / e L'^{flj) to be outside fij, we regard L^(rij) as a closed 
subspace of L'^{Q). The volume elements dV and dvj^'' of G and G^"' satisfy 
dV = Jjdvl°\ For A e B{L^{nj); L^{nj)), let A* and A*(j) denote their adjoint 
operators with respect to the volume element dV and dvj^^ , respectively. Then it 
is easy to show that 

A* = 

Taking A = Rf\z), and noting that R{z)* = R{z) and Rf\z)<^^ = Rf\z), 
we prove (jXTil) . By ((XTU)) and ((XT^ . EjB.jRf\z) and R{z)Vj{z)JjRf\z) are 
compact on which implies the last assertion of the lemma. □ 

3.2. Essential spectrum. 
Lemma 3.4. (Jess{H) = [0, oo). 

Proof. Lemma 3.3 implies XjR{z) — XjJj^Rf\z)JjXj is compact. Therefore 

N 

(3.15) R{z) = X,Ji'Rf\z)J,X, + K{z), 

where K{z) is a compact operator and satisfies 

(3.16) ||i^(z)|| <qimzr2(i + |^|), 

where |1 • || denotes the operator norm in L^(ri) and the constant C is independent 
of z. For /(A) G C^(K), there exists F{z) G C^(C), called an almost analytic 
extension of /, such that F{X) = /(A) for A G K and \d^F{z)\ < C„|Imz|", Vn > 0, 
and the following formula holds for any self-adjoint operator A : 

(3.17) fiA) = f d:F{z){z ~ Ay^dzdz. 

(See e.g. [25] or [30].) We replace (z - A}-'^ by ~R{z) and plug ([XTS]! . The 
inequality (|3.16[) implies \\dzF{z)K{z)\\ < C, and the integral over C converges in 
the operator norm, hence it gives a compact operator. We then see that (p{H) — 
J2f=iXjJi^f{Hf)JjXj is compact for any ip{X) G C^{R). Since cr(i?j"^) - 
[0, oo), we have ip{H^^^) = if ip{\) G C5"((— oo, 0)). Therefore ip{H) is compact if 
(p(A) G C^((-oo,0)), which implies that {-00,0) n cjessiH) = 0. For A G (0,oo) = 
cr(i/j°^), one can construct u„ G D{H^p) such that ||u„|| = 1, \\{H''p — A)u„|| 0, 
and suppun C {y > Rn} with i?„ 00. Then letting i;„ = Xj^n — SjBjUn, we 
have Vn G D(H), \\(H — A)f„|| — > 0, u„ — > weakly and ||wn|| > G uniformly in n 
with a constant C > 0. This implies A G (Tess{H). □ 
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The set of thresholds for H is defined by 

N 

(3.18) T{H) = \Jap{-AH,), 

where is the Laplace-Beltrami operator on Mj. Replacing Qq in §2 by ilj with 
j = 1, • • • ,N, we define the Besov type spaces Bj, B*. We put 

JV 

\\f\\B = \\xof\\mn)+J2\\Xjf\\B„ 

N 

MB'^\\xof\\LHn)+J2\\Xjf\\B;- 

The weighted space L"^'^ and the weighted Sobolev space iJ™'* are defined 
similarly. 

3.3. Radiation condition. A solution u e B* of the reduced wave equation 

( {H-X)u = f in O, A > 0, 

1 d,^u = on dCl, 
is said to satisfy the outgoing radiation condition if 



1 



where 



R 



(3.19) lim - / \\xj{dy - »P,(A))wi|i.(M^)dy = 0, 1 < < N., 







If dy — iPj{/^) is replaced by dy + iPj{X), we say that u satisfies the incoming 
radiation condition. In the following, u is always assumed to satisfy the boundary 
condition d^u = on dQ. 

Lemma 3.5. Let X G (0,cxd) \ T{H). If u & B* satisfies {H — X)u = and the 
outgoing (or incoming) radiation condition, it also satisfies 

1 

rToorJo \\Xj''\\lHM,)dy = 0, l<j<N. 

Proof. We take p{t) e Cg°((0,oo)) such that p{t) > 0, suppp(t) C (1,2) and 
/o°° p{t)dt = 1, and put 



V-R(y) = X (-|) , x{t) = / pis)d 



Then (pR{y) = 1 for y < i? and fR{y) = for j/ > 2R. We next construct 
V'H G Co°(f2) such that tf^R = 1 on /C and = on Q.j for I < j < N. Then we 
have 

u) = {i[H - A, u) = 0. 
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By the construction of 'i/'fl, [H,ijjR] — on IC. By the assumption (II. 2|) . on fij the 
commutator has the form 



(3.20) 



where Lj^r is a 1st order differential operator whose coefficients have the form 

^x(|)o(.-) 

and x{y) is either p{y) or p'{y). Let v — {1 + y)~^°u. Then by Lemma 2.1 and 
Lemma 2.3 (1) (which also holds for Agj), dxV,dyV G L^^^* for some < S < 1/2. 
Therefore 

I ^ {Wd.vWhiM,} + ll^y^^lli^CM,)) < {\\d.v\\ls + \\dyv\\ls) , 

which tends to as i? ^ cx). Therefore by Lemma 2.2 



lim (Lj,7?u,'u)i2(o.) = 0. 



Hence we have by using (|3.2Qp , 



(3.21) 



Assume that u satisfies the outgoing radiation condition. Using the inequality 



y 

R 
1 

R Jo 



1/2 



< C\\u\\b 
and (|3.2ip . we then have 

(3-22) }^^T.R p(|)(^.(A)x,":X,«)^.(M,)'^y = 0. 

j=i ° 

As in the proof of Theorem 2.7, we split XjU into two parts, 

XjU< = Ej{{-oo,X))xjU, XjU> ^ Ej{{X,oo))xjU, 

where Ej{-) is the spectral projection associated with — A/j^. Then by the short- 
range decay assumption of the metric. 

Since A ^ (t(— Am^ ), arguing in the same way as in the proof of p.7p . 

(-32 - A,,^. - X)xjUy = {-dl + S,)xju>, 

where is a self-adjoint operator on L'^{Mj) such that > 5{1 — Am^), (5 > 
being a constant. Therefore, Pj{X)xjUy G L^(rij), hence 
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Since Pj{X)xjU> — iCj{X)xjU>, where Cj{\) is a strictly positive operator on 
L'^(Mj), this also implies 

(3.24) ^rKi)ii^^">ii^^(-.)'^^-«- 

We show that 

In fact, in view of (|3.22p . splitting u = w< + u> and using (|3.23p . to prove (|3.25p 
we have only to show that 

and the same assertion with w< and u> exchanged. Let us note that 

< C'IIXj'«>||||Xj"<ll- 

Therefore 




since Xj'^< G ^* ■ By p.24p . this converges to 0. Similarly, we can prove p.26p 
with M< and u> exchanged. 

On the other hand, (Pj(A)xjU<, XjM<) > C'||xjU< ||^2(q .) for a constant C > 0, 
which depends on A. Therefore by (|3.25p 

(3.27) t^ooRj^ p(|)llx.«<lli2(Af,)rfy = o. 

By p.24p and (|3.27p . we complete the proof of the lemma. □ 

Lemma 3.6. Suppose u £ B* satisfies (H — X)u — f for A e (0, oo) \ T{H) and 
dyU = on dil. Assume also for some 1 < j < N , f G L'^''^{flj) for any s > 0, and 

(3.28) lim - / WxjuWl. ^M,)dy = 0. 

Then u £ L'^'^^^flj) for any s > 0. Moreover for any s > and any compact interval 
I C (0,oo) \ T{H), there exists a constant Cgj > such that 

(3.29) \\xM\L-..^n,) < Csji\\u\\B^n,) + ||/||l2.= (o,)), VA e /. 
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Proof. We construct counterparts of £j and Bj when the roles of G and cj"'' 

~(o) 

are interchanged. Namely, there exists an operator of extension £j such that for 
m > 1/2 and e H"^{dnj{R)) 

^ on dfljiR), 

on n\{njnn^{R-i/2)), 



supp (f/°V) c rij n fioo(-R - 1), 

f/*" e B(i7"'^*(af7j(i?));i?™+3/2'''(r2j-)), TO > 1/2, s > 0. 



If = on we have 

\j — y " "-"^i 'J 



a (o)(x,u) =2;-i-^»Bf^u on 



where 
We put 
Then 



c(0) _ ,,-l-eo?(0) 



Suppose u E B* satisfies {H — X)u = /, A e (0, oo) \ T{H), and d,yU = on dil. 
We put Vj — XjU — Sf'"' B^pu. Then satisfies 



(3.30) 



c?^(o)Wj=0 on (9i7j. 



Here Lj is a 2nd order differential operator with coefficients decaying like 0{y ^ ) . 
Note that /, := x,7 + L.^^u + (9^ + /^^^^)£f^ Bf \ e i2,i+.o(^.). 

Let Wj^n = {vj{-,y),ipn{'))L^iMj), where 'ipnix) is the normalized eigenvector as- 
sociated with the eigenvalue A„ of — A^^. Then we have 

(3.31) {-dy - Hn)Vj.n = gj,n, ^J.„ = X - A„, 

where gj.^ e L'^''^^+''°'>^^{{-oo,oo)), and Wj,„(?/) = 9j,niy) = for y < 0. Let 
ro(z) = (-9,2 - z)-i in ^^(R), i.e. 

/OO 
-oo 

where Imy^ > 0. Then as can be checked easily for any s > and 5 > 0, there 
exists a constant Cg^s > such that 

11(1 + |y|)Vo(-a)(l + |2/|)"'|1b(l2(k).l2(r)) < ^,5, Va > S. 

Therefore by (|3.3ip . one can show that 

(3.32) Ej{iX,oo))v, e L^^^^+'^^/^iflj), 
where Ej{-) is the spectral projection associated with — A/j^. 
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For A„ < A, we study Wj> separately. By (I3.28|) . 



R 



(3.33) lim - / \v,,„{yrdy = 0. 

In view of (j3.33|) . we see that satisfies both, of the outgoing and incoming 
radiation conditions. Adopting the outgoing radiation condition, we sec that vj^n 
is written as = ?'o(/in + iO)gj.n, i.e. 

vUy)^77^( f e^^(y~y"^g,M)dy' + r e^^^y'-y^gjAy')dy'). 

Note that 5j,„ e Li((0,oo)), since gj,„ G L2,(i+eo)/2((o, oo)). Therefore 

um v,.Ay) = \ e^^^y~y ^g,Ay')dy'- 

The condition (j3.33p imphes that this hmit is equal to 0. which implies 

v^y) - TT^ ( - / e^^^y-y'^g,,,{y')dy' + / e^^^y'-y^ g,^^{y')dy') . 

^yMn ^ Jy Jy ' 

Using the following Lemma 3.7 (Hardy's inequality), we have (1 + y)^^"^^-*^^?^ j,n G 
L2((0,oo)). Using we the have E L'^-''~^+'^''>^/'^{{Q,oo)). By Lemma 

3.2 and the formula XjU = vj + sf^sf^u, we have u E L2.(-i+«o)/2(f^^.)^ xhus 
we have seen that u gains the decay of order eo in Qj. Then in (j3.3ip . gj_n G 
L^'(i+^''')/2((0, oo). Hence Vj E i^'("^+^''/^)(%). Repeating this procedure, we 
obtain x^u E L^'™'°(rij), Vm > 0. The estimate p.29p can be proven by re- 
examining the above arguments. □ 

Lemma 3.7. Let f(y) E L"'^((0,oo)) and put 

/>oo 

uiy) = / f{t)dt. 
■'v 

The for s > 1/2 

y^i^-^)\u{y)\'dy<j^^-—^ y^V{y)?dy. 

For the proof, see [30] Chap. 3, Lemma 3.3. □ 

Lemma 3.8. Let arad{H) be the set of X ^ T{H) for which there exists a non-trivial 
solution u E B* of the equation [H — X)u = satisfying the radiation condition. 
Then (Jrad{H) — '^p{H) \ T{H). Moreover it is a discrete subset o/R \ T{H) with 
possible accumulation points in T{H) and oo. 

Proof. The first part of the lemma is proved by Lemmas 3.5 and 3.6. Let / be 
a compact interval in R \ T(H), and suppose there exists an infinite number of 
eigenvalues (counting multiplicities) in /. Let Un,n = 1,2,- •• be the associated 
orthonormal eigenvectors. 
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Take any R > and let xo — Xo be the function introduced in the beginning of 
this section. We decompose 

N 

Un = XoUn + ^(1 - Xa)X]Un- 

Then by p.29p . for any e > there exists R > such that ||(1 — XQ)un\\L^(n) < 
e uniformly in n. By the compactness of the imbedding of Hl^^{Vl) to L^^^(il), 
{XRUn]n is compact in L^(f7). Therefore {un}n contains a convergent subsequence, 
which is a contradiction. □ 

The set of exceptional points £{H) is now defined by 

(3.34) £{H)^r{H)[}ap{H). 

Weyl's formula for the asymptotic distribution of eigenvalues on compact manifolds 
and Lemma 3.8 imply that T{H) is discrete and £{H) has only finite number of 
accumulation points on any compact interval in R. 

3.4. Limiting absorption principle. For a self-adjoint operator H defined in a 
Hilbert space 7i, the limit 

\im{H - A =F ie)"^ =: (i? - A =F iO)"\ A G a{H), 

e—>0 

does not exist as a bounded operator on Ti. However if A is in the continuous 
spectrum of H, it is sometimes possible to guarantee the existence of the above limit 
in B(A'; A"*), where are Banach spaces such that X (Z H = H* (Z X* , and 

H is identified with its dual space via Riesz' theorem. This fact, called the limiting 
absorption principle, is central in studying the absolutely continuous spectrum, 
and many works are devoted to it. We employ in this paper the classical method 
of integration by parts pioneered by Eidus [16j . The crucial step is to establish 
a-priori estimates as in §2 of this paper, and to show the uniqueness of solutions 
to the reduced wave equation satisfying the radiation condition. After this hard 
analysis part, the remaining arguments are almost routine. 
We take any compact interval / C (0, oo) \ £{H) and let 

J = {z e C ; Re z e /, Im z 7^ 0}. 

We first note that Lemma 2.3 also holds for the solution to the equation 

r [H -z)u = f in 

I dyU — Q on 

by the standard elliptic regularity estimates. We put u = R{z)f and Vj = XjU — 
u as in the proof of Lemma 3.6, Then Vj satisfies (13.301) with A replaced by 
z. We can then apply Theorem 2.7 to see that 

(3.35) \\xM\b' <Cs{\\f\\B + \\u\\-s), 

for any s > 1/2, where C is independent of z e J. Once (|3.35p is proved, we can 
repeat the arguments in Chap 2, §2 of [3D] or those of Ikebe-Saito [57] without any 
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essential change. Note that here and in the sequel, we use ( , ) to denote the inner 
product 

{u, v) = uvdV 



Jn 

of L^(ri) as well as the coupling between B and B*, or L^'^ and L^'^'^. 

Lemma 3.9. Take s > 1/2 sufficiently close to 1/2. 

(1) There exists a constant C > such that 

snv\\R{z)f\\-s<C\\f\\B. 

(2) For any A G / and f E B, the strong limit linie^o ^^(-^ i ie)f exists in L^'^*'. 

(3) I 3 R{X ± iO)f e i^'"'' is continuous. 

Sketch of the proof. Suppose the uniform bound (1) is not true. Then there exist 
a sequence z„ £ J and /„ G B such that Zn = R{zn)fn satisfies ||u„||_s = 1 and 
^ 0- Without loss of generality, we can assume that 2„ ^ A G /. Using 
p.35p with < s' < s and the compactness of the embedding of Hf^^ into Lij^^, 
one can assume that u„ converges to some u G B* , and u satisfies the equation 
{H — X)u — and the radiation condition (see Corollary 2.6). Therefore w = by 
Lemma 3.8. However this contradicts ||u„||_s = 1. The assertions (2) and (3) are 
proved in a similar manner. □ 

Using this lemma one can prove the following theorem. 

Theorem 3.10. (1) For any X G /, liiHe^o -R(A ± ie)/ exists in the weak-* sense: 

3 lim(i?(A ± ie)f, g) =: (i?(A ± jO)/, 5), V/, geB. 
€—►0 

(2) There exists a constant C > such that 

\\R{X±iQ)f\\B.<C\\f\\B. AG/. 

Moreover R{X ± iO)f satisfies the outgoing radiation condition for A + iO and in- 
coming radiation condition for X — iO. 

(3) For any f,g £ B, I 3 X ^ {R{^ i *0)/i 5) continuous. 

(4) Let E(-) be the spectral decomposition of H. Then £'([0, 00) \ £{H))L'^{fl) 
= T-Cac{H), and we have the following orthogonal decomposition 

L'{n)^rLac{H)®'Hp{H). 

Sketch of the proof. Since L^'~^ [s > 1/2) is dense in B* , (1) follow from Lemma 
3.9 (2) and ((X^ . The assertion (2) follows from Lemma 3.9 (1) and ((X^ . The 
remaining assertions are proved in the same way as in Chap 2, §2 of 30j or Ikebe- 
Saito [H]. □ 

Let us recall that for a self-adjoint operator H = J°°^ XdE{X), the absolutely 
continuous subspace for H, HadFl), is the set of u such that {E{X)u, u) is absolutely 
continuous with respect to dX, and the point spectral subspace, 'Hp{H), is the 
closure of the linear hull of eigenvectors of H. 
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4. Forward problem 

4.1. Unperturbed spectral representations. Let {XjJjLo be the partition of 
unity defined in §3. Recall the spaces B and B* introduced in §2. For two functions 

f,gonQ, f — g means that 



lim — , 

-R— >oo it Jq 



" WxAv) ifi-^y) - 9{-.y)) \\l^M,)dy = o, i < Vj < n. 



We also use the same notation f ~ g for f,g defined on ^lj . 

Green's function of —Cp/dy^ — C on (0,oo) with Neumann boundary condition 
at y = is 

, I \ cos{^y)e'^y' , 0<y<y', 

VC [e'^y cosily'), 0<y'<y. 

Let Aj,i < Xj^2 < • • ■ be the eigenvalues of —^hj with normalized eigenvectors 
ipj,n{x)^ n — 1,2,- ■■. Without loss of generality, we assume that </5j,„ (cc) 's are 
real-valued. Let Hj^'' — ~dy — Af^. with Neumann boundary condition. Then 
R^p{z) — (Hj^^ — z)^^ is written as 



(4.1) [Rf{z)f){x,y) = Y. G{y,y'; 

n=l 



^ - ^J,n) iPj,nf) (x, y')dy' 



(Pj.nf) {X,y) = (/(•,2/),(/3j,„)^j>(x), 

( , ) being the inner product of L^{Mj ; det (hj) dx). Note that det(/iij) = 
detGf\ For }{x,y) e L'^{Mj x (0, oo); (det ''^dxdy) , we define its cosine 

transform by 

^cos (A)/(a;) =7r-i/2A-V4^" cos (yVx).f{x,y)dy. 
Lemma 4.1. For f E B, and A G (0,oo) \ (Tp(— A/i^. ), we have 



Proof. We first show that the right-hand side of (|4.ip is a bounded operator from 
B to B* . The sum over the terms in which Aj.„ > A is rewritten as 
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where fj,n = Pj.nf and kn = sj ^j,n — A. Then we have 

\\M^)f{;y)\\hiM) 

poo 

<Cx \{f{-,y),Vj,n)\^dy 



Aj,„>A ' 
|2 



< C'a 11/11 i2(MjX(0,Oo)) • 

Hence Aj{X) e B(L2; l2(-^^.^-) ^ b(S;S*). To estimate the term in which 

^j,n < ^, we put 

poo 

'>J'j,n{x) = / G{y, y'; A ± iO - Xj,n)fj,n{^^ y')dy' . 
Jo 

Then we have 

/•OO 

\uj,n{x)\<Cx \fj^„{x,y)\dy. 
Jo 

Since 

We have proven that Rf\\ ± iO) e B{B;B*). 

Now the assertion of the lemma is easy to prove if there exists no > such that 
fj^n = for n > no, and fj^n is compactly supported for n < no- Since such an / 
is dense in S, we have proven the lemma. □ 

The generalized eigenfunction of J?j°^ is defined for A > Xj^„ 
(4.2) 'i>f,l{x,y;X) = 7t-'/\X - Xj,nr'^^ cos (y,/X - A,>)<p,- 
This ^j^l^{x, y; A) is often denoted by ^f^iX) in the sequel. It satisfies 



(4.3) 



(-A^(o, -A)vl'(°^(A)=0 in Qj, 

5.(o)*;5(A) = o on anj. 



The Fourier transformation associated with H^^^ is defined by 
(4.4) (A)/ = £c,^,„(A)^j2(A)/, 

n=l 

where Cj,„(A) is the characteristic function of the interval (Aj,„,oo), and 



(4.5) 
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where dvj°^ = {detCf^Y'^cLxdy. Define a subspace of {{0, oo); L^{Mj);dX) by 

OO 

(4-6) ri . 

= w2jnW^j,n{x) ; /„ e L^((Aj,„,oo) ; c'A) I . 

Tlien defined by = ^(A)/ for / e Co°°(%) is uniquely extended 

to a unitary operator 

We put 

N 

(4.7) h = 0L2(Af,), 



wfiere L^(Mj) = L'^{Mj; ^ydeijhj) dx) , and and also 

(4.8) ^(°)=.(^r,...,^(,°)). 

By the computation similar to the one to be given in the proof of Lemma 4.3 
below, one can show that 

^ ([i?f (A + zO)-i?f (A-^0)] /,/) = ll^f (A)/||i.(M^.). 

Therefore, T^°\x) G B{B;L^{Mj)), and J^j"^(A)* G B(L2(Mj); K*). 

Here we must pay attention to the following remarks. The first one is that in 
(1131), T^/\X) is a finite sum: 

(4.9) ^f(A)= J2 

The second remark is that the adjoint * is taken in the following sense: 



(4.10) (.Ff (A)/,Mi.(M,) = (/,^f (A)*/i)i.(n,)= / f:F^°\xrhdVj'\ 

Jfij 

(h e L^iMj )). Therefore 

(4.11) .Ff (A)*= ^j'nW*^ 

Aj,„<A 

and for h G L^{Mj) 

(4.12) {:F^"lW*h) {x,y) = M>\^l{x,y;X)iK^,,„)LHM,y 
Since J^j°'(A)* satisfies (ff"^ - A)J^]°^(A)* = 0, we have 

TfiXy eB{L^M,);H'^-n, s > 1/2, 

hence 

TI°\X) G B(ff-2-^;L2(Af^.))^ 5 > 1/2. 
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4.2. Perturbed spectral representations. Using Bj and Vj{z) in Subsection 
3.1, for A > Xj n we define the generalized eigenfunction for H by 

(4.13) ^,,„,±(A) = ix, - £,B,) *;°1(A) - R{\ T iQ)Vj{mfl{^)- 

Here putting s = (1 + eo)/2, we regard £jBj and V^j(A) in B{H'^-^]L'^''). Note 
that \l/j,n,±(A) G B* . This definition easily implies 



(4.14) 



(-Ag-A)«',-„,±(A) =0 in fl, 
9^*j,n,±(A) = on on. 



The generalized Fourier transformation for H is defined by perturbing JFj"-* . We 
put for A > Aj,„ 

(4.15) ^,,„,±(A) - 4°!(A)J,(x, - {E.B^y - V,{Xr)R{\ ± zO), 

where J, = (det G/det Note that {£jBj)\Vj{X)* e B(L2'-''; i/-2^'*), and 

R{\ ± iO) e B(L2.^; ij2.-s) p B(i/-2,^. x2.-«), hence (jiJSll is well-defined. 

Lemma 4.2. For / G C5"(ri), we /laue 

(4.16) (.F,,„,±(A)/)(x) - (^J^¥~^fdV^ ^,,,,{x), 

1 /2 

where dV = (det (G)) dxdy. 

Proof. We put w = (xj - (^jSj)* - V,(A)*)i?(A ± iO)f. Then by using (1410)) 

(•?^j,n,±(A)/,/l)L2(M,) = {^jjlWJjU,h)L2(^M,) 



uTfl{\YhJ,dvl'\ 
We then use (|4.12p to see that the right-hand side is equal to 



= ((X, - {£,B,r ~ y,(A)*)i?(A± zO)/,vI/(°)(A)) ()i:^^.(M,) 

= (/,*j,n,±(A)) {ipj.n,h)L2(^M,), 

which proves the lemma. □ 
The adjoint operator ^j,n.±(A)* is defined by the following formula: 

(4.17) (.^,.„,±(A)/»L2(M,) = {f,J'j,n.±Wh*)mn), h e ^'(M,). 
Lemma 4.3. The adjoint operator J-'j^„,±{X)* has the following expression: 

(4.18) ^,,„,±(A)* = (x, - - i?(A T *0)T/,(A)) .Ff2(A)*, 
where the adjoint J-'j'^(A)* is taken in the sense of J^. 
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Proof. Let u = {xj - {£jBj)* - Vj{X)*)R{\ ± iO)f. Then as is shown in the 
proof of Lemma 4.2, 



(^,,„,± (A)/, h)L. (M, )^ in ^j°„) (A)* h J,dvj°^ 



Plugging the form of u, we see that the right-hand side is equal to 
(/, (x, - £,B, - i?(A T iO)V,{X))T^'^{Xrh)L.^n), 

which proves the lemma. □ 
We define 

oo 

(4.19) ^,;±(A) = ^c,-„(A).F,-„,±(A)= ^ ^,,„,±(A), 

n=l Aj,„<A 

(4.20) ^±(A) = (.Fi,±(A),--- ,^Ar^±(A)). 
Lemma 4.4. For any X G (0, oo) \ £{H) and f £ B, we have on flj 



(4.21) R{X±iO)f~±iV^ - A,- „)-i/'*e±W^-^-".Fj,„,±(A)/. 

Proof. This follows from (|3.14p . Lemma 4.1 and the definition (|4.15p . □ 

Lemma 4.5. For any X E (0, oo) \ £{H) and f E B, we have 

^ mx + zO) - R{X - zO))/, /) = \\T±iX)f\\l 

Proof. We prove the case for J^_|_(A). We have only to prove the lemma when 
/ G C^{fl). We compute in a way similar to that in Lemma 3.5. Take p{t) G 
C^((0, oo)) such that p{t)dt = 1, and put x(i) = p{s)ds. Let u = R{X+iO)f 
and 

AT 

V'i? = xo + Yx{-^)xj{y), 

where {xj}jLo is the partition of imity on fl, and y in Xjiv) is the local coordinate 
on flj. We then have 

{[H - A, iIjh]u, u) = {ipRU, /) - (/, TpRu). 
As M G yB* , by computing the commutator [iJ, ■f/^i^] , we then have 

rI^J}!^ {p{-^)Xjiy)dyU,uj = {u,f)-(f,u). 
i=i 

Since u = R{X + iO)f satisfies the radiation condition (see Theorem 3.10 (2)), 
{dy — iPj{X))xjU ~ 0. Therefore 

N 
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Now we note that 



Let v± be the term in the right-hand side of (|4.2ip . Using Lemma 4.4, we first 
replace u of the right-hand side of {Pj{X)u, u)^2(^m ) ^^^^ move Pj{X) to 

the right-hand side of the inner product, and replace uhy v±. Since Pj n{X)ip^^^^ — 
y/X— Aj,„ ^1^1, we have Pj(A)Jv,-,„^+(A) = ^ A — Aj,„ J^j,n,+ (A) . The lemma then 
follows from a direct computation. □ 

The formula in Lemma 4.5, when integrated with respect to A over (0,oo), is 
a counterpart of the Parseval formula in the Fourier transformation, and a crucial 
step for the spectral representation. Using Tlj in (|4.6p . we put 

N 

(4.22) n^^n,. 

The following theorem can be proved in the same way as in [28^ or Chap. 3 of [30] . 
Theorem 4.6. (1) For A ^ T{H), JP±(A) G B{B;h). 

(2) The operator (.F±/)(A) = ^±{X)f defined for f G B is uniquely extended to a 
partial isometry with initial set TiadH) and final set Ti.. 

(3) (.F±iI/)(A) - A(.F±/)(A), VAe (0,(^)\f(i/), yfeDiH). 

(4) ■^±(A)* G B(h;S*) is an eigenoperator of H with eigenvalue A in the sense 
that 

{H - X)T±{X)*i' = 0, VV'Gh. 

(5) For any compact interval I C (0, oo) \ T{H) and g G H, we have 



J'±{X)*g{X) dX G L^{n). 

Let In be a finite union of compact intervals in (0,c») \ £{H) such that /„ C In+i, 
U$^i-^n = (0,oo) \ £{H). Then for any f G TiaciH), the inversion formula holds: 

/ = s-lim / T±{XY{T±f){X)dX. 



turbed Laplacian in the end ilj . 



4.3. Time-dependent scattering theory. Let Hj = —dy — /S.h. be the unper- 



Theorem 4.7. The wave operator W± : 0,=ii^(^) ~^ L^{fl) defined by 



N 



W+=s-\imye''^Xje~''^"^ 

i— »±oo — ^ 



exists and is complete, i.e. RanTU± — HadH). Moreover 
(4.23) W± = {T±yT'-°\ 

where T'^^^ is the Fourier transformation defined by 14-S^ for the system of Lapla- 
cians {h[°\--- ,H^°^). 
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Sketch of the proof. We argue in the same way as in Chap. 2, Theorem 8.9 of 
[3D]. Take / e HadH) such that „,+/) (A) e C^f ((A,- „, oo)) and .Fj,„,+/ = 
except for a finite number of n. Then by Theorem 4.6 and Lemma 4.3 



e 



= / e-"^.F+(A)* {j^+f) {\)d\ 





oo 



(4.24) ,.00 



^/ (^j°i(A) (^,,„.+/)(A)dA 

e-^'^%B,[T^%)) (.F,,„,+/)(A)dA 



J./m 

pOO 



J2 / e-*^i?(A-zO)y,(A) (.Fj,':^(A))* (.F,;„,+/)(A)(iA. 

j,n 



Because of the decay of £j, the 2nd term of the right-hand side tends to in L^(Jl). 
Letting A — ^/7l, we have 

{H-k^ + iO)-^ = {A-k + iO)-\A + k)-\ 

We then put 

g{k) = 2kiA + k)-%ik') (^^^(fc^))* (.F„„,+/)(fc^), 

g{t) = / e-**'=g(fc)dfc. 
Jo 

We show that 

(4.25) \\git)\\ <Cil+t)-^~\ t>0. 

In fact, take h E L'^{fl) and consider 



{g{t),h)^ 2ke-^''^{ (Tfl{k^)) {T,,r.,+fW).V,{k^){A + kr^h)dk 
= 51 / ^^/"^ |^e^«(*fe+aV'='-^^.") + g-«(tfe-aV'='-^j-.")j ...dk. 

Since V,(fc^) contains a factor (1 +?;)~^~^, sphtting the integral suitably and inte- 
grating by parts, we have 

\{g{t)M<C{l + t)-^-^\\h\\, 

which proves (I4.25p . We use the notation f{t) ^ g{t) if \\f{t) — g{t)\\ — > as i ^ 00. 
In view of the following Lemma 4.8, we obtain as t — > 00 
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+ • 



in L^(ri). This implies the existence of the hmit 

N 

(4.26) s - hm y e'V^f x.e""^ = f J"^"^) * ^- 

Since (jF^*')) JF_(_ is a partial isometry with initial set Ti.ac{H) and final set L^(il), 
(|4.26|1 also implies for g = (51, • • • ,gN) e ®f^iL^{^j) 



AT 



(4.27) 



Let us prove this fact. We put U{t) = J^f^i e"^"' Xje~"^- Then ffTI^ 
implies that U{t) (•^") ^+ —'■ U strongly, which implies 



(4.28) 

We show that 



U{t)* U* weakly. 



N 



(4.29) ||C/(0*5ll-||.9ll = l|C/*5ll, 5 = (ffi, • • • , 5a) G 0i'(%). 
In fact, we have 

A 



A 



A 



^l|2 II -itJH'-"' I 



By the scattering property of e , || (1 — Xi)e 

A A 

i=i j=i 

Now, (|4:28l) and ([QS]) yield ||C/(0*5 - U*g\\ 0. This completes the proof of 
Theorem 4.7 for W+. The assertion for is proved similarly. □ 

Lemma 4.8. Let A be a self-adjoint operator on a Hilbert space Ti. For f{k) e 
Co((0, 00); 7i), we put 



gjW — > 0, which proves 



f±{t) 



±ikt 



f{k)dk. 



Then for any e > 



{A-kTiey^e^'^*f{k)dk 



< 



\f±{s)\\ds. 



Proof. This is proved in [5D], Chap. 2, Lemma 8.10. For the reader's convenience, 
we reproduce the proof. By virtue of the identity 



/•OO 

{A-kT *e)^^ e'^'^^^-'^^'^Us, 
Jo 



we have 



/ {A-kTier^e'^''''f{k)dk = ±i e^''^^^'''> f±{s + t)ds, 
Jo Jo 

which proves the lemma. 



□ 
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4.4. S-matrix. The scattering operator is defined by 5* = W-. We consider 

its Fourier transform : S = T'^"'> S 

Lemma 4.9. We have a direct integral representation: 

{Sf){x)^s{x)f{x), yfeii, VA>o, 

where S{X) = (^Sjk{X)) -^^^ ^^^^ is a bounded operator on h called the S-matrix, and 
is written as follows 

SjkiX) = 5jk - 2^iT,^+{X)Vk{X) (t^\x) 

Proof. Lemma 4.5 implies 

(i?(A + iO) - R(X - iO)) = J'±(X)*T±(X). 

2m 

By Lemma 4.3, we tlien have 

Tk,+ {XY - Tk,-W* = 2TiiJ^+{XrT+{X)Vk{X)Tf\x)*. 
Then we have by Theorem 4.6 (2), for f,gEH 

{{T+-T^){T+rf,g)^^2mY, / [nX),T+{X)V,{X)(T'^\X)) g{X)) dX. 

By S ^T+{T-y . Hence the lemma follows. □ 

Let hj(A) be the linear subspace of L'^{Mj) spanned by ipj^n such that Aj,„ < A 
and put 

N 

h(A) = 0h,(A). 

Then 5(A) is a partial isometry on h with initial and final set h(A). The scattering 
amplitude is defined by 



Let Aj^^kniX) : L^iMk) ^ i^(Mj) be given by 

(4.30) A,™,fe„(A) - J^,.^,+ {X)Vk{X) [4°lWy ■ 
We then have 

Sjk{X) ~ 5-jkIj = ~2m ^ Ajyn,kn{X), 

where Ij is the identity operator on L'^{Mj). When j, k and the energy A > is 
fixed, (Ajm,fen(A)) is a finite matrix of size {dj,dk), where dj — 4/^{{j,rn); Xj^m < A}. 
Let Ajm,kn{X) be defined by 

(4.31) Ajm,kn{X) = {Ajm.kn{X)ipk,mfj,m)i^2(^]yj^-j ■ 

Then we have 

(4.32) Ajm,kniX)h = -4jm,fen(A)(/l, Vk,n) L^(Mk)Vl,ra, V/l G L'^{Mk). 
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The scattering amplitude is computed from the asymptotic expansion of the 
generahzed eigenfunction in the foUowing way. 

Lemma 4.10. 



p 



(A-A,,„)V4 

Proof. This directly follows from (|4.13p and Lemma 4.4. □ 

5. From scattering data to boundary data 

5.1. Non-physical scattering amplitude. In this section, we observe waves 
coming in from and going out of the end fJi assuming that 

(5.1) Gi = {dyf + hi{x,dx) on Vli. 

This amounts to studying the scattering amplitude Aim^i„(A) of (|4.30p . which is 
rewritten as 

Al,n,ln{X)=T,os{\-Kn.)Pl.nJl{xi-Vl{\yR{\ + lQ)) 
■Vi{\){T,os{\-Xi.n))* Pl.n- 

Note that Bi — Q, because of the assumption (|5.ip . By the expression (|3.12p . 
Vi(A) and Vi(A)* are independent of A and compactly supported in the ?/- variable. 
Therefore, ^im.in(A) defined for A > max {Ai Ai_„} is analytically continued to 
the upper half plane C+ — {Ini A > 0}. This analytic continuation can be extended 
to a continuous function on C+ U (M \ 8{H)). We denote the obtained function for 
{A < max{Ai,„, Ai,„}} \ £{H) by ^^'^'^^^(A) and caU it the non-physical scattering 
amplitude. These functions can be represented by (|5.2p . where Tcos{^ — M,m) and 
•?^cos(A — Ai.n) are replaced by their analytic continuations. Let 

(5.3) <^fl{x,y-\) = ^-i/2e-^^/4(Ai,„ - \)-^'^cos\, (y^Xi,n - a)(^i,„(x), 

and put, similarly to (|4.5p 



(^cosh(Ai,„ - A)Pi,„/) [x) = $S(A)/dF/"^^ ViA^). 

In the following, we always assume that X^£{H). The explicit form of A^^i|^(A) 
is given by the following lemma. Recall that the non-physical scattering ampli- 
tude A']^'1^^{X) coincides with the physical scattering amplitude ^im,i„(A) for 
A > max {Ai,,n, Ai_n}. 

Lemma 5.1. (1) If Xi^m < X < Ai,„, 

ATinW = ^cosiX - Ai,„,)Pi,„Ji(xi - ViiXyRiX + tO)) 

■VliX) {,Tcosh{Xl,n-X))* Pi,n. 

(2) IfXi^n < A < Xi^„i, 

^wl(A) = T,osh{Xi^,n - A)Pi,™ Ji(xi - ^i(A)*i?(A + zO)) 
•l^i(A) (.Feo.(A-Ai,„))*Pi,„. 
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(3) IfX< niin{Ai,™, Ai,„}, 

= ^cosh{Xi,m - X)Pi.mJi{xi - ViiXyRiX + iO)) 
■Vi{X) (^cosh(Ai,„-A))*Pi,„. 

In accordance with (I4.13P , we define non-physical eigenfunction by 
(5.4) $i,™,±(A) = xi^^'LW - T *0)yi(A)<i>(°UA). 

Note that the physical eigenfunction ^i^m,-(A) defined for A > Xi^m is analyti- 
cally continued through the upper half space C+ to the nonphysical eigenfunction 
*i'i.m,-(A) defined for A < Ai.m- The non-physical scattering amplitude is computed 
from the asymptotic behavior of non-physical eigenfunction in the following way. 

We put 

Then we have for h e L'^{Mi) 

Lemma 5.2. (1) //Ai,„i < A < Ai,,i, we have as y oo, 
(2) If X < max{Ai.m, Ai^„}, we have as y ^ oo, 

Pl.m (A) - $i;,(A)j ^ _ ^^1/4 l„(A)V'i,„, 

with a super exponentially decreasing error, that is, with the error r(y) satisfying 
\riy)\ < CNC-^y for any N > 0. 

Proof. The assertion (1) is proved in the same way as in Lemma 4.8. By (|4.ip . 
letting C = A — Ai^m, we have as y ^ oo 

ie'-^y r°° / ^ \ 

Pl,™i^l(A-^^0)/(a;,y)~ cos [./Cy'j Pi,,nf{x,y')dy' 

with a super exponentially decaying error. This, together with p.l4p and Lemma 

5.1, proves (2). □ 

5.2. Splitting the manifold. We take a compact hypersurface F C i^i having the 
following property. 

(C-1) F splits fl into a union: Q = Clext U flint so that flext H flmt = F, flint is a 
manifold with smooth boundary, and Qext C fli- (See figure 2.) 

Let O C flint be an open, relatively compact set such that it has a smooth 
boundary not intersecting dflint and that flint \ O is connected. Denote flQ — 
flint \ O and 

r F if = 0, 

^'^ ~\dO if O 7^ 0. 
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Figure 2. Surface T splits fl to two parts, manifold flint with a 
smooth boundary and its complement Clext C r^i. 

We put for f,ge L^iTo) 

dSx being the measure induced from the metric G on Vq- We put Ho = — in 
r^o endowed with the Neumann boundary condition: 

(5.5) dyV — Q on ^^lo^ 

V being the unit normal to the boundary. If has only one end, ^int is a bounded 
region. If f2 has more than one end, is unbounded and the spectral theory 
developed for H applies also to Hq. To see this, we have only to replace K, by 
/CU(rii n ^int \ and to argue in the same way as in §3 and §4. Let £{Ti.o) be ap{Ha) 
when f2o is bounded, and the set of exceptional points when Qo is unbounded. 
Next we consider the case O = so that To — T. 

Lemma 5.3. Suppose A ^ £{H) U S{H^), and let \l/i,„__(A) and <i>i,n,-(A) be 
physical and non-physical eigenf unctions for H . Then the linear subspace spanned 
by 9,y\E'i,„,-(A)|p, 9y<i>i_„^_(A)|p, n = 1, 2, . . . , is dense in L'^{T). 

Proof. We show that, if / G -^^(r) satisfies 

(5.6) (/,a,*i,„,_(A))r = (/,9.$i,„_(A))r =0, V?i > 1, 

then / = 0. We define an operator S'^ e B{{H^^^{T)y;H-^{n)), where {H^^^{T)y 
is the dual space of if ^/^(F), by 

and put u = R{X — iO)S^f by duality. This means that , if G_ (A; X, X') is Green's 
function, i.e. the integral (Schwartz) kernel of i?(A — iO), 

u{X) = {R{X-iO)S'^f) [X)^ j^d,,G-{\-X,X')f{X')dSx'. 

where d^i means the conormal differentiation with respect to the variable X' . Then 
u Cz B* , and by (j3.14p . we have the following asymptotic expansion on fJi 

a.(A)e-'V^-^i." (/,a,*i,„_(A))j,¥.i,„(x). 

Ai,„<A 
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In particular, if Xi,n < A 

(5.7) (u, (^i,„) ~ C„(A)e-'^\/^^ (/, a,M'i,„,_(A))p , 
C„(A) being a constant. In a similar way, we have for Ai_„ > A 

(5.8) (u, (^i,„) ~ C;(A)e-V^i."-^ a, (A))j, 

modulo a super exponentially decaying term. Note that m„ — (u, <pi,„) satisfies the 
equation {—dy + Ai_„ — A)m„ = for y > a, a being a sufficiently large constant. 
In view of the assumption of (|5.6p and (|5.7p . (|5.8p . we then have (u, (pi,„) = for 
y > a, hence u(x,y) = for y > a. The unique continuation theorem then implies 
u = on flext- By the property of classical double layer potential, d,jU is continuous 
across F, so that d^u\r = 0. 

Next we show that u = in Clint- In the region flint, we have (—Ac — A)u = 0. 
If Qint is bounded, then u = since A is not a Neumann eigenvalue. If Qint is not 
bounded, u satisfies the incoming radiation condition, since so does u in Q. Then 
w = in flint by Lemma 3.4. As u — R{X — iO)(5p/ G Lf^^{il), it follows from the 
above that w = in il. Applying iJ — A, we have (5p/ = as a distribution, hence 
/ = on r. □ 

5.3. Interior boundary value problem. For z E C \ 8{Ho), we consider the 
following boundary value problem 

{{Ho — z)u = Q in r^o, 
d^u^Q on dVLo\To, 
d,u = f eHl'^{To) on To- 

The incoming radiation condition is also imposed, if flint is unbounded and z G K.. 
The Neumann-Dirichlet map (N-D map) is then defined by 

(5.10) Ko{z)f = u\^^, 

where u is the solution to (|5.9p . When O = 0, we use for the N-D map of the 
operator i/g the notation Kq{z) = A(z). 

Now we consider the operator theoretical meaning of the N-D map. Note that 
from now on O may be a non-empty set. We put T = (J-"c, J-p), where J^c is the 
generalized Fourier transform for Ho (which is absent when flint is bounded) and 
Tp is defined by 

Tp :Hp{Ho) 3ui-^ {{u,'ipi),{u,ip2),- ■ ■), 

and where 'Hp{Ho) is the point spectral subspace for Hq and ipi is the eigenfunction 
associated with the eigenvalue \i of Hq ■ There are two kinds of generalized Fourier 
transformation, J--\- and J--. Both choices will do as Tc- Then is a unitary 

(5.11) T ■.L''{fi,nt)-^n®<c\ 



MANIFOLDS WITH CYLINDRICAL ENDS 33 

where d = d\m'Hp{Ho)- li d = oo, = . Moreover, we have 

(5.2) ,«„-„-..^"wr^,,,|,_A_, 

where are the eigenprojections associated with eigenvalues A^, numbered count- 
ing multiplicities by i = 1, 2, . . . , d, and the right-hand side converges in the sense 
of strong limit in L'^iflo). 

Let rro G B{H^{no); H^/'^{To)) be the trace operator to To, 

We define Sto € B((i7i/2(r^))'; as the adjoint of rr^: 

{Srof.w)mn,^,) = {f,rro^)LHro)^ f ^ (i?'/'(ro))', w G H\no). 
With this in mind we write 

Lemma 5.4. For z ^ £{Ho), the N-D map has the following representation 

J* ^,(A)*.F,(A)<5r„ ^ ^ J^S*RSra 



Jo A - z 



= 1 ^ 

Proof. For / e Hy^iTo), take / e W^iVLo) such that a^,/ = / on To and 
/ has compact support in fimf. Then the solution u of (|5.9|) is written as u = 
/ - (ffo - z)-1(-Ag - z)/. Let .9 = .F,(A)(Ag + z)/. Then for any h e H,„t, 
where Hint is defined by (|4?22l) with j = 2, . . . , iV. 

(.Fc(A)(Ag + z)J, h) = ((Ag + z)f,Tc{\rh) 

= idJ,rraJ'c{Xrh)LHro) + ili^G + z)TAXrh) 
= (/,rro.Fe(A)*/i)i2(r^) + (/, {-X + z)TAXr h). 

This implies 

.Fc(A)(Ag + = .F,(A)(5ro/ + (-A + z)T,{X)l 

Hence 

.F,(A)*.Fe(A)(AG + z)/^^ ^ .Fe(A)*.F,(A)Jro/ ^, _ 



Similarly, 



X — z In X — z 



^ P.(Ag + z)/ _ ' P.Srof _ ^ p 7 

2—1 i—1 i—1 



Since J^^J^c/ + Eti -P^/ = by (lET^ . these imply that 

which proves the lemma. □ 
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Let US call the set 

(5.13) [s*r^Tc{XrTc{X)Sra-A & iO,^)\£iHo)} U [[K,6*^^P,Sro)Y^^^, 

where d = dim TYp (iJg) ) , the boundary spectral projection (BSP) for Ho on To- 
On the other hand, the set 

(5.14) {Tc{X)Sto ; a G (0, oo) \ £{Ho)} U { (a„ V,(x) | j,^) }^^^ 

is called the boundary spectral data (BSD) on To- 

By using the formula (|3.17[) . we have the following lemma. 

Lemma 5.5. For a bounded Borel function ip{X) with support in M\T{Ho), where 
T{Ho) is defined by 113. 18\) with j = 2, • • • , N , we have 

S*r^^{Ho)Sra ^ / ^{X)S*r„T,{XrT,{\)SrodX + Y,viX,W^oP^^^o■ 

Proof. By the formulae (|3.17p and (|5.12p . this lemma holds for any (p(X) E C^{M.\ 
T{Ho))- The general case the follows from the approximation. □ 

Usually BSD is referred as given data in the inverse boundary value problems. 
What is actually used in our reconstruction for the manifold is the BSP. 

Lemma 5.6. Let O C ^int- Then knowing the N-D map Ao{z) for all z ^ ij{Ho) 
is equivalent to knowing the BSP for Ho - 

Proof. By Lemma [5^ one can compute the N-D map by using BSP. Taking tp{X) 
as the characteristic function of the interval [a, t) and taking note of the remark 
after p.34p . we differentiate the formula in Lemma 5.5 with respect to t to recover 
(^f o-^c(t)*^c(t)'5ro for t e M \ £{Ho)- Since 

p 5l^J^,(XYTe{X)5ro 
= ^o{z) - / ^— dA, 



one can obtain eigenvalues A^ as the poles of the right-hand side. The residues 
in these poles provide us with (Jpg X^a =A^j'^rc>- This determines the terms 
^Y^Pj^To for indexes j such that Xj = Xi, up to an orthogonal transformation 
of the eigenspace associated to the eigenvalue Ai, see [311 Lem. 4.9] or [33]. Thus 
we can determine the BSP for Ho- □ 



We complete this section by the following result used later to prove Theorem ll.il 
Let OW, r = 1,2, be as in Theorem ll.il We take F as above, which moreover has 
the following property: G^^^ = G^p on Vlext — ^^ext — ^e^t- We put the superscript 
(''^ for all relevant operators and functions explained above. Let A('")(A), r = 1,2 
be the N-D map for h''^\ that is, when = 0. The basic idea of the following 
Lemma is due to Eidus [IB] . 
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Lemma 5.7. Under the assumptions of Theorem ll.ll we have A(^'(A) = A'-^'(A) 
for A e (0,oo) \ Ur^i^2{£{H'^'-'>) U £{hI^'^)), and BSP's for H^^^ and H^^^ coincide 

on r. 

Proof. Since ^^^•'(A) = s[^i{X), the physical scattering amphtudes coincide, 
hence so do non-physical scattering amplitudes by analytic continuation. Let 
u = *1^1,_(A)-*1^^ _(A) andw = $S,_(A)-$J^,^,^_(A). Then since H^i) i/^^) ^ 
—dy — Ah^ on Qext, u and v satisfy (— 9^ — Ah^ — X)u — and {—dy — Af^ — X)v = 
in Qext- Using Lemma 5.2 and arguing in the same way as in the proof of Lemma 

(r) (r) 

5.3, we have u = u = in flext- Therefore, '^i „ _ and ^\ „ _ as well as their 
normal derivatives coincide for r — 1,2 and for all n £ Z+ . Since they satisfy the 
equation ((5Jl) for i/0 = H^''\ we have A'^^^X) = A(^)(A) due to Lemma 5.3. The 
last statement now follows immediately from Lemma 15.61 □ 

6. Boundary control method for manifolds with asymptotically 

cylindrical ends 

In this section we reconstruct the isometry type of the manifold {Qmt, G) using 
given data. 




Figure 3. We will construct manifold flint by iterating local con- 
structions. First, a neighborhood Ui C flint of S C F is recon- 
structed. Next, a ball O = B{Xi,p) C Ui is removed from the 
manifold and data analogous to measurements on dO are con- 
structed. After that, the metric is reconstructed in a larger ball 
B{Xi,t), and the procedure is iterated to reconstruct the whole 
manifold flint- 



Theorem 6.1. Assume that we are given the set T as a differentiable manifold, the 
metric G on T , and the ESP for Hii,. These data determine the manifold {flint, G) 
up to an isometry. 

For proving this theorem, we use the boundary control method described for 
compact manifolds in e.g. [7l[32l[35]. The reconstruction of non-compact manifolds 
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is considered previously in the conference proceedings 34J and in 8 with different 
kind of data, using iterated time reversal for solutions of the wave equation. 

The proof of Theorem 16. II is divided into a series of lemmas. Our reconstruction 
of {^lint,G) is of recurrent nature. We will begin with the case when O = so 
that we are given just the set To = F as a differentiable manifold, the metric on it, 
and the BSP for the operator iJg on T. We apply the boundary control method to 
reconstruct the metric G on some neighborhood Ui of F. Then, we will take a point 
Xi e U\r and p > such that B{Xi,2p) C C/i, where B{Xi,r) denotes the ball 
of radius r with center at Xi. We take O = B{Xi,p) and show that we can find 
the BSP for the operator Hq on F^ = dO. Then we apply the boundary control 
method starting from Tq , which would allow us to recover {Qint , G) in a larger 
neighborhood U2 D Ui of F. Proceeding in this way, we will eventually recover 
the whole of (r2int,G). Therefore, our further considerations deal with arbitrary 
O C flint including the case = 0. 



6.1. Blagovestchenskii's identity. Let us first consider the initial boundary value 
problem 



Lemma 6.2. Assume that we are given the set To as a differentiable manifold, 
the metric G on Tq and the BSP for Hq on To ■ Then for any given f,h£ 
G^{To X M-|_) and t, s > these data uniquely determine 



(6.1) 





Proof. Let 



S{t,X) 



sin(\/Ai) 



Then the solution uf{t) is written as 
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Using the similar decomposition for u''(s) and the fact that J-c{fJ,)J-c{X)* — (5(// — A), 
we obtain the foUowing formula: 

(6.2) 

pt pS n OO 

= / dt' I ds' I d\S{t-t\s-s\\){5-^^T,[\rT,{\)5rof{t'),h{s'))^,^^^^ 

d 



"'0 "'0 



/ dt' [ ds'J2sit-t\s-s',X,){S*r^P,drof{t'),h{s')) 
Jo Jo 



where S{t, s, X) — S{t, X)S{s, A). Observe that the right-hand side depends only on 
BSP and the metric on To- □ 

Above, the formula (|6.2p is a generalization of Blagovestchenskii identity (see 
[32( Theorem 3.7]) for non-compact manifolds. 

6.2. Finite propagation property of waves. Let us next introduce some nota- 
tions. For t > and E C To arbitrary, let 

r!o(s,t) ^{x erio; do{x,^) < t} 

be the domain of influence of S at time t. Here, do{X,Y) is the distance between 
X and Y in fio. We use also the notation na{Y,t) = flo({Y},t). More generally, 
when / = {{'Sj,tj)}'j^-^ is a finite collection of pairs (Sj, ij), where C and 
tj > 0, we denote 

,7 

^o{I) = U ^o{^j,tj) ^{X eno ; do{X, E,) < for some j = 1, . . . , J}, 

For any measurable set B C ilo, we denote L'^{B) — {v ^ L'^{no); w|no\B = 0}j 
identifying functions and their zero continuations. 

Lemma 6.3. Assume that we are given the set Vq as a differ entiahle manifold, the 
metric on Tq and the BSP for Hq on Tq. Then, for any given f € C^{Tq x R+), 
T > 0, and I = {(Sj, <j)}^^j^, where Yij C Vq o,fs open sets or single points, and 
tj < T, we can determine 

(6.3) aj,Tif)^ f \uf{T)\'dV. 

Jno\no{i) 

Proof. When S C Fq is an open set and h g C^(S x ]R_|_), it follows from the 
finite velocity of wave propagation (see e.g. Sec. 4.2], see also [3(L Ch. 6]) that 
the wave u^{t) = w''(- ,t) is supported in the domain ilci(E,t) at time t > 0. It 
follows from Tataru's seminal unique continuation result, see [HIHH], that the set 

(6.4) {u^\t); h e C^(S X R+)} 
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is dense in L^(rici(S, t)), see e.g. [32l Theorem 3.10]. This clearly implies that, 
when T > and / = {(Sj, where Sj are open and tj < T, the set 

Xf := {u^'iT); h = hi + --- + hj, G C^{^j x [T - tj,T])} 
- span,=i,...,j {u^t,); h G Cq^CE, x [0, t,])} 

is dense in i^(rio(/)). 

Next, we consider the non-linear functional 

aiM = inf{||^i^'-''(r)||i.(,,„); h = h, + ■ ■ ■ + hj, h, e C^iY., x [T - t„T])} 

where / € C^(ro x R+), T > 0, and / = {(S^-, C To are open, and 

tj < T. By the formula (|6.2p . the BSP and the metric on To determine the value 
a/_r(/) for any /. Moreover, as uf^^{T) = {T) — u^(T) and Xf is dense in 
L'^{Qo{I)), we see that 

(6.5) ai,T{f) = 11(1 - Xnoii)W{T)\\Uno)' 

where xno (/) (^) the characteristic function of the set flo (I) on ■ This proves 
the lemma for the case when all Ej are open. 

(k) 

If for some j, the set Ej is just a point Xj e To, we define for those j's E^- C 
To, fc = 1,2,--- to be open neighborhoods of Xj such that E^'^'*'^-' C E^'^'' and 
Pi J, E^'^' = {Xj}. For those j's for which Ej is open, we define E^'''' = Ej. Denote 
the corresponding finite collection of (E^'^^tj) by /(fc). Then 

oo 

no{i{k + 1)) c no{i{k)), noil) = f] no{i{k)), 

fe=i 

and for any b € L'^{ilo), 

(l - Xoo(/(fe)))^ ^ (1 - Xf2o(/))^, a.e. as fc ^ oo. 

As 1(1 - Xno(i(k}))b{-)\ < 1(1 - Xnoii))K-)\^ a.e., using the monotone convergence 
theorem, we see that 

ai{k).T{f) ^ II (l - Xno{i)W {T)\\\-,(no) = «/,r(/). 
Thus, the BSP and the metric on To determine aj^rif) for such I's. □ 

Definition 6.4. Let / = {(Ej, /' = {(E^-, t^)}/^i and T > 0, where 

Ej,E^- C To and tj,ij- < T. We say that the relation / > /' is valid on manifold 
flo if 

(6.6) rJo(/') \ n,o{I) has measure zero. 

Lemma 6.5. Let I = {(Ej, tj)}/^i, /' = {(E^ , and T > 0, where Ej, E^- C 

To are open sets or single points and tj,tj < T. Assume that we are given the set 
To as a dijferentiable manifold, the metric on To , the BSP for Ho on To , and the 
collections I and I' . Then we can determine whether the relation I > I' is valid on 
manifold Qo or not. 
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Proof. The relation / > /' is valid on manifold flo if and only if 

(6.7) aiM<ai',Tif) ior idl f e (To x R+). 

Indeed, the equivalence of (|6.6p and (|6.7p follows from (|6.5p and the fact that, by 
Tataru's density result (|6.4|) . the functions u^{T), f £ C^{ro x K+), are dense in 
L'^(no{To,T)). As for given /, by Lemma [6T3l we can evaluate both sides of (|6.7p . 
using the BSP and the metric on To, these data determine, for any pair {1,1'), if 
the relation / > /' is valid or not. □ 

For any Xq G ^lint \ dflint , introduce the exponential map 

expxo : 7(Xo, 

where f G S'xo(a„0 = {?7 G Tx„in,nt);\v\ = 1} and < t < s(Xo,C). Here 
7(Xo,C) (^) is geodesic on ft, parametrized by the arclength, with 7(Xo,4) (0) = ^o, 
7(jC(,^^)(0) = and [0, s(Xo,^)) is the maximal interval of t, when 7(Xo,{)(i) stays 
in flint, that is, s(Xo,0 = sup{i ; 7(Xo,C) ([0, i)) C fi^t \ af^„t}. Denote by 

(6.8) s(Xo)= inf s(Xo,0 
so that 

5(Xo,5(Xo)) ca„t\5r! 

int • 

Define now 



TiXo,0 = sup {t; d0(7(Xo,c)(<),^o) = t}. 

0<t<s(Xo) 



At last, define 



(6.9) r(Xo)=. inf r(Xo,C). 

In geometric terms, the above definition of t(A"o) means that in the ball B{Xo, t{Xo)) C 
f^int \ dQint , it is possible to introduce the Ricmannian normal coordinates 

X^{^,t):^e Sxo {n,„t),0 < t < t{Xo) 

which satisfy 7(Xo.C)('') — 

We also need the boundary exponential map 

expr„ : {{Z,t) e To x R+ ; < i < so(Z)} 3 {Z,t) ^ 7(z,.)(<) e flo- 

Here v is the interior unit normal (with respect to flo) to Tq and 

(6.10) so{Z) = sup{i > ; 7(z,.) ((0, t)) C flo \ dflo). 
For any Z € To , let 

(6.11) To{Z)= sup {i;do(7(2,,)(i),ro) = 0- 

0<t<so(Z) 

In the following, we impose the following condition (C-2) on S. 

(C-2) For O = , E is an open subset of F such that d^{Y., dV) > 0, and for O ^ 0, 
S = do. 
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We define 

(6.12) To{^) = ugToiZ). 

In geometric terms, the above definition of re)(S) means that, in the set 

i(S, ro(S)) = {7(z,.) (i) ; Z e S, < i < ro(S)} c {flo \ dflo) U S, 
it is possible to introduce the boundary normal coordinates 

X ^ {z,t),z Gi:,o<t< Tc.(s) 

satisfying X = j(z,u){t). Observe that when O = B{X,p), X e flint \ flint and 
p > is small enough, then 

ToidO) = r{X) - p. 

Lemma 6.6. Assume that T, c To satisfies condition (C-2). Let Y G T,, Z G 
To, t < re)(S), and X = j(Y,v){t)- Assume that we are given the set Fq as a 
differentiable manifold, the metric on To and the BSP for Ho on To ■ Then we 
can determine the distance do {X, Z) on flo ■ 

Proof. Note that as t < tc)(S), the set rio(y, t) \ no{To,t — s) contains a non- 
empty open set for all e > 0. For s,e > 0, let us denote (see Fig. 4) 

I,{t) = {{Y,t),{To,t-e)}, lUt,s) = {{Z,s),{To,t-e)}. 



X 




Y T Z 



Figure 4. In the figure O = $ and s is small enough so that the 

sot Q,o(Y,t) \ floiToTt — e) is not contained in Q,o{Z,s). This is 
the situation when I'^{t, s) ^ Itit). 

Let us next show that for any r > there is eq > such that 

fto{Y,t)\no{To,t- e) C B{X,r), when e < eq. 

If this is not true, there are r > 0, a sequence Sj 0, and Xj e floiY,t) \ 
flo^Xoi t — ^j) such that do{Xj, X) > r. As flo{Y, t) is compact, by considering a 
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subsequence, we can assume that Xj converge to X G rici(y, t). Then 

do{XX) = lim do{Xj,Y)<t, 

do{X,To) = Hm do(X,,To) > t, 

implying that y is a closest point of Tq to X and c?o(X, Y) — t. Let us recall that 
the shortest curve from a point in to Tq, which end point is an interior point 
of To, is a normal geodesic to Tq- Thus, we see that X = 7(y,,y)(t) = X, which is 
in contradiction to d{X,X) > r. Thus, the existence of Sq for any r is proven. 

The above implies that when s > d{X,Z), the set ilo{Y,t) \ floiToit ^ £) is 
contained in Q,o{Z, s) for all sufficiently small e > and therefore, 

(6.13) there is ei > such that /^(i, s) > /^(t) for all < £ < £i. 

On the other hand, for s < d{X, Z), the set flo{Y,t) \ f^o(ro, t - e) 7^ do not 
intersect with no{Z,s) at any e > small enough and thus (|6.13p does not hold. 
Thus, by Lemma 16.51 we can find doiX, Z) for any Z g as the infimum of all 
s > for which (^1^ holds. □ 

For S C Fo satisfying (C-2) and < T < rci(E), let iVs,T and A/s,t be the sets 

Afs,T = {Xer!o; X = 7(y,.)(t), 0<t<T, y eS}, 

(6.14) 

Ms,T = {Xeno: X = 7(y,.)(<), < < < T, y e S} C iVs.T = Ms,T. 

Note that AI-^^t is open in fla. 

6.3. Boundary distance functions and reconstruction of topology. Let us 

next consider the collection of the boundary distance functions associated with 
To- For each X e flo, the corresponding restricted boundary distance function, 
rx G C{Vo) (note that Vq is compact) is given by 

rjf:Fo^M+, rx{Z) ^ do{X, Z), ZeVo- 

The restricted boundary distance functions define the boundary distance map Ro ■ 
flo C{To), Ro{X) ~ rx- The boundary distance representation of N^^t C ilo 
is the set 

i?o(A^E,T) = {rx e C(Fo); X e N^,t}, 
that is, the image of iVs,T in Ro- Clearly Rr^ '- ^ Ci^o) is continuous. 

Lemma 6.7. Assume that we are given the set To as a differ entiable manifold, the 
metric on Tq, the ESP for Hq on Tq, an open set S C Tq satisfying condition 
(C-2), and < T < ToiTi). Then we can determine the set 

Ro{N^.t) = Ro{{i(Y,.){t); yes, 0<t< T}). 

Proof. By Lemma |6^ for y e S, t < T and Z g Tq, we can find do{X, Z) where 
X — 'y(Y,v){'t) from BSP. This gives us the function rx{Z), Z G Tq, and for such 
X's- Thus, BSP and the metric on To determine the set Ro{M^^t)- Using (|6.14p . 
we obtain i?o(A^s,T) by closure of Ro{M^^t) in CiVo)- D 



42 



HIROSHI ISOZAKI, YAROSLAV KURYLEV, AND MATTI LASSAS 



Consider properties of Ro- Assume that rx — ry for some X,Y E A^s.t- Let 
Z G Tq be the point where the function rx attains its minimum. Then, it is the 
closest point of Tq to X. Thus, the shortest geodesic from X to Z is normal to 
To, i.e. X — "l{z,v){i) with t — rx{Z). The same arguments show that Z is also 
the closest point of To to Y and t — ry(Z), and hence Y — 7(z.i/)(0- Thus X = Y 
and Ro is injective on A^s,t- 

Thus, map Ro '■ -/Vs.t ^ Ro{Ny:,t) is a bijective continuous map defined on 
a compact set, implying that it is a homeomorphism. This implies that the map 
Ro ■ Ms^T Ro{M^^t) is a homeomorphism. As BSP and the metric on To 
determine the manifold Ro{Ms,t) with its topological structure inherited from 
C{To), we see that these data determine the manifold M^^t as a topological space. 

Lemma 6.8. The set Ro{M^^t) C C{To) can be endowed, in a constructive way, 
with a differentiable structure and a metric tensor G, so that {Ro{M^^t)iG) he- 
comes a manifold which is isometric to {M-^.t, G) with Ro being an isometry. 

For compact manifolds, the result analogous to Lemma 16.81 is presented in detail 
in [32l Sect. 3.8]. Since the proof is based on local constructions, it works for 
non-compact manifolds without any change. However, for the convenience of the 
reader, we present this construction. 

Proof. Let us define the evaluation functions, Ez, Z e To, 

Ez : i?o(Ms,T) ^ K, Ez{rx) = rx{Z) = do{X, Z). 

For r(-) e RoiMs,T) corresponding to a point X S M^^t, i-e. r(-) = rxi-), we 
can choose points Zi, . . . ,Zn G To close to the nearest point of To to X so that 
X {do{X, Zj))"^^ forms a system of coordinates on fto near X, see [321 Lem. 
2.14]. Similarly, the functions Ezj, j = l,...,n, form a system of coordinates 
in Ro{M^^t) near rx- These coordinates provide for Ro{M^^t) a differential 
structure which makes it diffeomorphic to manifold My,^t- 

Let us denote by G the metric on Ro{M-s,t) which makes it isometric to {M^^t, G), 
that is, G = {{Roy^)*G. Let r e i?o(Afs,T) and X e Ms,t be such that r = rx- 
Let Zq is a point where r obtains its minimum, that is, the closest point of To to X. 
When Z is close to Zq, the differentials of functions Ez are covectors of length 1 on 
{Ro{M^,t), G), see [351 Lem. 2.15]. This is equivalent to the fact that the gradients 
of the distance functions X ^ do{X, Z) have length one. By this observation, it 
is possible to find infinitely many covectors dEz, Z G To of length 1 at any point 
r of Ro{My.,t)- Using such vectors, one can reconstruct the metric tensor G at 
r. By the above considerations, BSP determines the manifold {M^^t,G) up to an 
isometry. □ 

6.4. Continuation of the data. Let us now consider the case when = and 
we are given the set F as a differentiable manifold, the metric G on F, and the BSP 
for . Assume that there are two manifolds fj'^j and Sl^-^j such that F is isometric 
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to subsets r^J) C for j = 1,2 and that the BSP for iJ^^', j = 1,2, comcides 

with the given data. Let now S C F satisfy condition (C-2) and 

0<r<min(r«(S),r(^)(E)). 

Then the above constructions show that the manifolds 

M^^^ = {Xe X = 7(y,.)(i), < i < T, y e E} 

with j ~ I and j ~ 2, are isometric. Thus, we can consider the set M^}p, denoted by 
Ui as a subset of both manifolds fi^^j and fi-^j, and, by the previous considerations, 
we can construct a metric G on it which makes {Ui,G) isometric to (M^lp, G^^'^), 
J = l,2. 

We continue the construction by continuation of the data using Green's functions, 
cf. [Mllini- To this end, let z e C\R+ and consider the Schwartz kernel Goiz; Y, Y') 
of the operator {Hq — z)~^. It satisfies the equation 

(6.15) {Ho-z)Go{z;-,Y')^5y', YX ^ ^ ^.nt\0, 

d,Go{zr X)\di-io -0. 
We denote G(z; Y, Y') = Go{z : Y, Y') when O 0. 

Lemma 6.9. Let U C fii„t he a connected neighborhood of an open set E C F, 
where E satisfies condition (C-2) with O — $. Let Xq (z U \ dflint and p > be 
such that O — B{Xq,p) (1 U \ d^lint- Assume that we are given the metric tensor 
G in U . Then BSP on F for the operator iJ^ determines G{z; Y, Y') for Y,Y' Cz U 
and z G C \ £{H(d). Moreover, these data determine BSP on Tq for the operator 
Ho- 

Proof. By Lemma [5^ BSP on F determines the N-D map A(z) at FxF. By Lemma 
15.41 the Schwartz kernel of the N-D map A(z) at F x F coincides with G{z] Y, Y'). 
Thus we know the function G{z;Y,Y') for Y,Y' e E. As the Neumann boundary 
values of y i-^ G{z;Y,Y') on F \ {Y'} vanish, using the Unique Continuation 
Principle for the elliptic equation (|6.15p in the Y variable, we see that the values 
of G{z;Y,Y') are uniquely determined for F' G E and Y € U \ {Y'}. Using the 
symmetry G{z;Y,Y') = G{'z;Y' ,Y) and again the Unique Continuation Principle, 
now in the Y' variable, we can determine the values of G{z;Y,Y') in {{Y,Y') G 
U X U; Y Y'}. Considering G{z;Y,Y') as a locally integrable function, we see 
that it is defined a.e. in U x U. 

For Y' G (noOU) \ dno, denote by G'g\z;Y,Y') a smooth extension of 
Go{z\Y,Y') into O. Then 

(-Ag - z)G'g\z- Y, Y') - 5{Y, Y') = F{Y, Y') G G^{n„,t), 

where suppF(-,y) C O. Therefore, 

Go{z;Y,Y') = G{z;Y,Y')+ f G{z;Y,Y'')F{Y" ,Y')dVY" . 

Jo 
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In particular, 

(6.16) 9,(r)G(z; Y, Y') + / a,(y)G(z; Y, Y")F{Y'\ Y')dVY" =0, Ye dO, 

Jo 

where viY) is the unit normal to O at Y. On the other hand, if F{- ,Y') e C°°{U), 
suppi^(- , y ) C O, satisfies (|6.16p . the function 

(6.17) Giz;Y,Y')+ [ G{z;Y,Y")F{Y",Y')dVY", Y,Y'eU\0, 

Jo 

is Go(z; Y, Y'). As we have in our disposal G{z; Y, Y') for Y, Y' E U, we can verify 
for a given F, condition (|6.16p . 

Now, we return to ^i -^j, fi -^^^ with F and BSP on F being the same. We denote the 
associated functions appearing above by adding the superscript (j). Let (|6.16p holds 
with G{z:Y,Y'), F{Y",Y') replaced by G^^^z;Y,Y'), F^^^Y",Y'), respectively. 
Since G^^\z;Y,Y') = G'^^^z;Y,Y') onU xU, KWi also holds with G{z;Y,Y'), 
F{Y",Y') replaced by G'^^'> {z:Y,Y'), F^^^Y",Y'), respectively. Thus, for Y,Y' e 
U\0, we have 

g'^1\z;Y,Y') = G^JHz;Y,Y') + [ G'^^\z-Y,Y'')F'^^\Y"X)dVY", J = 1,2, 

Jo 

so that 

G^^^(z; y, Y') = Gg'(z; Y,Y'), z e C \ R, y, y' G C/ \ O. 
In particular, this implies that A^-'(z) — K^q{z), 2; e C \ M. Then by Lemma [5?6l 
BSP's for and h'^^ coincide. □ 

Next we show that we can use these data to determine the critical distance which 
we use in the step-by-step construction of the manifold. 

Lemma 6.10. Let Xq e \ dQint and < p < t{Xo)/2. Let O = B{Xq,p) 
and To — dO. Assume that we are given the set To as a differ entiahle manifold, 
the metric G|p^ on To, and the BSP for Ho on To- Then these data determine 
To(ro) = -r(Xo) - p. 

Proof. Let us assume that t^ < t(A'o) — p. Then, for any Y e Fo, the set 
^^0(^,^0) \ ^o{^0,to — £) contains an open neighborhood of 7(y,^)(io — £/2) and, 
therefore, has positive measure. Hence, \it< t{Xq) — p, then the condition 

(6.18) Vy e Fo Ve > : := {{To,t-e)} ^ L'y^^ {{Y,t)} 
is valid. 

Let us next assume that condition (j6.18p is valid and consider its consequences. 
First, observe that by (|6.8[) and (|6.9[) . we have either 

(a) s{Xq) = t(A'o) and there is y G Fo such that X = ^(y.v){t{Xo) — p) G dilint, 
or 

(b) s{Xq) > t{Xo) and there are Y e To and s such that s{Xq) > s > t{Xq) — p 
and do{'-f{Y,u){s),To) < s. 

Let us consider these two cases separately. 
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(a) It follows from ()6.8|) and ()6.9|) that X is a closest point to Xq on dilmt- 
Therefore, the geodesic 7(f.i/) intersects dflint normally a.t X = 7(y,i/)(s), s = 
riXo) ~ p. 

Assume next that t > is such that 

(6.19) Ve > : It,e t 4,f 
Then for any e > there is 

(6.20) e f7o(F,i)\rio(ro,t-e). 

As r2o(Y,t) is relatively compact, there are £„ ^ and X„ = X^^ such that 
Xn X' ^ ^int as n — > GO. Then 

(6.21) doiX',Y) = t, do{X',To) = t. 

This shows that Y is the closest point of Tq to X' in ■ Consider a shortest curve 
fi(s) from Y to A"'. By [2], a shortest curve between on a manifold with boundary 
is a C^-curve. Moreover, it is a geodesic on il.a\dil.o- Since ^{s) is a shortest curve 
from X' to dflo, it is normal to dil,a at Y. Thus /i(s) = 7y,iy(s), s < t(Ao) — p. 
However, j(Y.u)is) hits dflint normally at s = r(Ao) — p. Therefore, by the short-cut 
arguments, we see that the curve 7(y_i^)([0, r(Ao) — p]) C flo can not be extended 
to a longer curve which is a shortest curve between Y and its other end point. Thus 
p C 7(y,,^)([0, t(Ao) — p]), implying that t = do{Y,X') < t{Xo) — p. Hence in the 
case (a) the condition (|6.18p implies that t < t{Xo) — p. 

(b) In this case arguments are similar but slightly simpler. Again, assume that 
t > is such that (|6.19p is satisfied. Again, there are e„ > and A„ = X^^ 
satisfying (|6.20p . such that A„ — > X' and X' E flint satisfies (j6.2ip . Moreover, a 
shortest curve p{s) from Y to X' coincides with the normal geodesic 7(y,i/)(s) for 
small values of s. Since the geodesic 7(f,j/)([0, s']) is a shortest curve between its 
end points for s' < t(Ao) — p but not for s(Ao) — p > s' > t{Xq) — p, we see that 
M C 7(y,i.)([0, t(Ao) - p]) and thus t < t{Xo) - p. 

Therefore, in both cases (a) and (b), the condition (|6.18p implies that t < t(Ao) — 
p. Combining these facts, we see that 

t{Xo) — P = sup{t > ; condition (|6.18p is satisfied for t}. 

The lemma then follows from this and Lemma 16.51 □ 

6.5. Proof of Theorem 16. 1[ We are now in a position to complete the proof of 
Theorem 16.11 

6.5.1. Local reconstruction of Riemannian structure. We start our considerations 
with = 0. Let E C F satisfies condition (C-2) and T > be a sufficiently small. 
In fact, we can consider any < T < T0(S). Using Lemma l67fl we see that the set 
i?0(Ms,T) C C(r) is uniquely determined. On this set we introduce the boundary 
normal coordinates, 

r{-)^{Z,t), t= mm r{Z'), 
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where Z is the unique point on E on which r(-) attains its minimum. Observe that 
these coordinates on Ri^^M-^ x) coincide with the boundary normal coordinates of 
the point X G ^int such that 

r{-)=rx{-). 

Thus, Rit,{MY.,T) with the above coordinates is difFeomorphic to My.,t- 

Next we use Lemma l6.8l to endow R(i,{My.^t) with Riemannian metric, G, so that 
{R(Ii{M^^t)iG) is isometric to manifolds {M^^TiG). 

Remark. For the inverse scattering problem considered in the introduction, Sec- 
tion 6.5.1 is not necessary, because we know a priori the Riemannian structure of 
the open set (f2i„t \ d^mt) H Sli. However, to make the results of §6 appropriate 
for general non-compact manifolds with asymptotically cylindrical ends, we have 
included this step. 

6.5.2. Iteration of local reconstruction. To describe the procedure which we will 
iterate, let us assume that Ui C Hint is a connected neighborhood S C F which 
satisfies condition (C-2) with = and that we know the Riemannian manifold 
{Ui,G) up to an isometry. Since the set {R(d{Mj:^t,G)) is already determined, we 
can take Ui = M-^^t, where T > is sufficiently small . 

Choose Xi G Ui and p > such that O = B{Xi,p) C C/i. By Lemma [6.91 we 
can determine G^(y, Y') for all Y, F' e C/i and z e C \ M. Then Lemma [Ql gives 
us BSP on do. Therefore by Lemma [6.10i these data determine tc)(Fo), hence 
t{Xi) = ro(Fo) + p. Take any X e B{Xi,t) \ O, where r = t(Xi), and let Y 
be the intersection of dO and the geodesic with end points Xi and X. Taking any 
Z G do and applying Lemma |6.6[ we can then find do{X, Z). 

Using, similarly to the above, Lemmas 16 . 71 and 16 .81 we can find the image of the 
embedding Rq ■ B{Xi,t) \0 ^ C{dO). We then recover, in the boundary normal 
coordinates associated with dO, i.e. the Riemannian normal coordinates centered 
at Xi, the metric tensor G on B{Xi,t) \ B{Xi,p), and, since G on B{X,p) is 
known, on the whole B{X,t). This construction makes it possible to introduce 
the structure of the differentiable manifold on C/i |J B{Xi, r) which we considered, 
by now, as a disjoint union of two Riemannian manifolds. Next we glue these two 
components together. To this end we observe that, since O C Ui, we have in our 
disposal Green's function G{z; Y, Y') for F, F' e O and z e C\R. The set O can be 
considered also as the subset B{Xi,p) of B(Xi,t), and thus we know the function 
G(z; Y, Y') for Y, Y' g B{Xi,p) e.g. in the Riemannian normal coordinates centered 
at Xi. Thus, using the Unique Continuation Principle, we can determine, in the 
Riemannian normal coordinates, the function G(z; Y, Y') for all Y e B{Xi, r) and 
Y'eB{Xi,p). 

Since Y' i— > G{z; Y, Y') is a smooth function in flmtXiY} and G{z; Y, Y') ^ oo as 
y Y, we see that for Yi, 1^2 S ^int, we have Yi = Y2 if and only if Gz{Yi,Y') = 
GziY2,Y') for all Y' e Hint, z e C\M. Using the Unique Continuation Principle, 
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this is equivalent to G{z;Yi,Y') = G{z;Y2,Y') for all Y' e B{Xi,p), z e C\R. 
Next, let us define that the points Xjj G C/i and Xg E B{Xi,t) are equivalent 
and denote Xu Xg if G{z;Xu,Y') = G{z;Xb,Y') for all Y' e B{Xi,p), z G 
C \ K. Then the manifold U2 ^ UiU B{Xi, r) C flint is diffeomorphic to manifold 
{Ui U B{Xi,t))/ ^, which is obtained by glueing together the equivalent points on 
Ui and B{Xi,t). As we know the metric tensor on both Ui and B{Xi,p), we have 
reconstructed a Riemannian manifold (C/2, G) C {^int, G) up to an isometry. 

6.5.3. Maximal reconstruction. Let us iterate the above process, that is, we start 
from an open set S C F satisfying condition (C-2) with = 0, construct its 
neighborhood Ui, and iterate the construction by choosing at each step j — 1,2, . . . 
a point Xj G Uj and constructing a Riemannian manifold isometric to J7j+i = 

UjUB{Xj,t{x,) cn,nf 

Consider the open sets in ^int \ d^int which can be reconstructed, with the 
metric, when we are given the set F with its metric and the BSP on F. As the 
collection of these sets is closed with respect to taking the union, consider maximal 
open set Umax C flint \ dflint which can be reconstructed, with its metric, from the 
set F with its metric and the BSP on F. Let us show that Umax — flint \ dflint- 
Since flint \ dflint is connected, it suffices to show that Umax is open and closed. 
By construction, Umax is open. Let now X ^ dflint be a limit point of Umax, i-e., 
X = lim„^oo Xn, Xn G Umax- Denote a = d{X, dflint) so that if ^ G a/4), 
then s{Y) > 3a/ A, see (|6.8p . Since the cut locus distance of the Riemannian normal 
coordinates is continuous with respect to the center, see e.g. [Ml Sec. 2.1] or [3T], 
there is J > such that t{Y) > S for all Y e B{X, a/4). 

Let now Xn G Umax satisfy the inequality d{Xn,X) < a = min(a/4, (S/4). 
Let us assume that Xn has a neighborhood B{Xn, Pn), with a sufficiently small 
Pn < d{Xn,X), which can be reconstructed using N(n) iteration steps, that is, 
B{Xn,Pn) C UN(n)- Then T{Xn) > 4(T so that X G i3(X„,r(X„)). By Lemma 
16.91 we can find the BSP for the operator Hq with O = i?(X„,r(X„)) and, us- 
ing one more iteration step, to reconstruct the Riemannian structure on UjsK^n) U 
B{Xn,T{Xn)) which includes the point X. Therefore, the point X is in Umax- 
This shows that Umax is relatively open and closed in flint \ dflint- Thus, Umax — 

fliiit \ dflin- 

The above shows that using an enumerable number of iteration steps we can 
construct a Riemannian manifold isometric to (flint \ dflint, G). Thus we have 
reconstructed the Riemannian manifold {flint \ dflint, G) up to an isometry. 

It remains to identify the differentiable and Riemannian structures near dflint - 
Observe that flint is just the closure of flint\d flint with respect to the distance func- 
tion generated by the metric G on flint \ dflint- Moreover, for any open relatively 
compact set E c dflint, 7'0(S) > S > 0. 

Let < t < S and consider the set 



St = {X G fl^nt \ dflint d{X , dfl,nt) = t, d{X , Z) = t, ioi SOmC Z G S}. 
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This implies that for X G St the closest point Z G ilmf is in S and X — ^z,v(i\ 
Therefore, St is a smooth {n — l)-dimensional open submanifold in r2,;„t of points 
having the form X = J{z,v){t), Z E J^. This makes it possible to introduce the 
boundary normal coordinates in Ms, a which provides the differentiable structure 
near E. Writing the metric tensor G in these coordinates and extending this ten- 
sor continuously on E, we find the metric tensor in ilint in the boundary normal 
coordinates associated to S. □ 

6.6. Proof of Theorem 11.11 Having Theorem 16.11 in our disposal, Theorem 11.11 
follows immediately from Lemma 15.71 □ 
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